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ABSTRACT 


The feasibility of using satellite-derived thermal data to generate 
realistic synoptic-scale winds within the planetary boundary layer (PBL) is 
examined. Diagnostic ’'modified Ekman” wind equations from the Air Force Global 
Weather Central (AFGWC) Boundary Layer Model are used to compute winds at seven 
levels within the PBL transition layer (50 m to 1600 m AGL) . Satellite- 
derived winds based on 62 predawn (0921 GMT 19 April 1979) TIROS-N soundings 
are compared to similarly-derived wind fields based on 39 AVE-SESAME II rawin- 
sonde (RAOB) soundings taken 2 h later. Actual wind fields are also used as a 
basis for comparison. Qualitative and statistical comparisons show that the 
Ekman winds from both sources are in very close agreement, with an average 
vector correlation coefficient of 0.815. Best results are obtained at 300 m 
AGL. Satellite winds tend to be slightly weaker than their RAOB counterparts 
and exhibit a greater degree of cross-isobaric flow. The modified Ekman winds 
show a significant improvement over geostrophic values at levels nearest 
the surface. 

Horizontal moisture divergence, moisture advection, velocity divergence 
and relative vorticity are computed at 300 m AGL using satellite-derived winds 
and moisture data. Results show excellent agreement with corresponding RAOB- 
derived values. Areas of horizontal moisture convergence, velocity convergence, 
and positive vorticity are nearly coincident and align in regions which later 
develop intense convection. Vertical motion at 1600 m AGL is computed using 
stepwise integration of the satellite winds through the PBL. Values and patterns 
are similar to those obtained using the RAOB-derived winds. Regions of maximum 
upward motion correspond with areas of greatest moisture convergence and the 
convection that later develops. 

ii 


ACKNOWLEDGMENTS 


The authors gratefully acknowledge the many helpful 
suggestions of Drs. Gandikota V. Rao> Yeong— jer Lin and David 
O'C. Starr. Special thanks go to Mr. Gary Jedlovec for his 
assistance in data acquisition; the Technical Services Branch 
of AFGWC for providing information on the Boundary Layer Model; 
the USAF Environmental Technical Applications Center Library 
for background research; and particularly to Mr. Mark Fenbers 
for assisting in preliminary calculations, typing the manuscript 
and drafting the figures. 

This research was sponsored in part by the National 
Aeronautics and Space Administration under Contract NAS8-33370 
and is under the auspices of the Atmospheric Sciences Division, 
Systems Dynamics Laboratory, NASA Marshall Space Flight Center, 
Alabama. The author is also indebted to the Air Force Institute 
of Technology for providing a scholarship to Saint Louis 


University . 


TABLE OF CONTENTS 


Title Page 

Acknowledgements Ill 

Table of Contents Iv 

List of Tables vll 

List of Figures vlii 

1. INTRODUCTION 1 

a. Background 1 

b. Statement of the problem 5 

2. REVIEW OF THE PBL THEORY 8 

a. PBL structure 8 

b. Classic surface layer formulations 9 

c. Classic Ekman layer formulations 13 

d. Higher-order closure models 16 

3. THE AFGWC BOUNDARY LAYER MODEL 19 

a. Fundamental aspects of the model 19 

b. Surface layer equations 20 

c. Computational scheme for the surface layer 23 

d. Transition layer equations 26 

e. Computational scheme for the transition layer . » . . 28 

4. DATA 31 

a. Satellite data 31 

b. Radiosonde data 34 

c. Synoptic situation 36 

5. METHODOLOGY 40 

^ ^ Gcnsir 2 ti.on of i.npu't d3.t3s 

iv 


40 


TABLE OF CONTENTS (CONTINUED) 


Page 

b. Objective analyses of the Input data 40 

1) Grid structure 40 

2) Temperature and dewpoint 42 

3) Pressure 42 

c . Computation of modified Ekman winds 44 

d. Computation of standard geostrophlc winds 44 

e. Validation of technique 45 

f . Response determination 45 

g. Computation of kinematic parameters 48 

6. RESULTS 51 

a. Comparison of Input data 51 

b. Evaluation of RAOB-derlved winds 61 

1) "Actual" winds 61 

2) RAOB-derlved winds 68 

3) Validation 88 

c . Satellite-derived modified Ekman winds 88 

d. Sensitivity analysis 100 

e. Kinematic parameters 112 

7. SUMMARY AND CONCLUSIONS 125 

APPENDIX A - DERIVATION OF MODEL EQUATIONS 129 

APPENDIX B - RESPONSE DETERMINATION 146 

APPENDIX C - CORRELATION COEFFICIENTS 150 

APPENDIX D - SCALE ANALYSIS OF HORIZONTAL MOMENTUM 

EQUATIONS 151 


V 


TABLE OF CONTENTS (CONTINUED) 


Page 

APPENDIX E - BAROCLINICITY ^ THE PBL 153 

Bibliography 157 


vi 


LIST OF TABLES 


Table Page 

1 Average and maximum values* as well as standard 

deviations (a) , of the gridded RAOB and satellite 
temperature, dewpoint and pressure fields 57 

2 Mean values and standard deviations (m s~^) for 

"actual" wind data 72 

3 Same as Table 2, except for RAOB-derived standard 

geostrophlc and modified Ekman winds 81 

4 Standard deviations of "actual" u, v and wind 

speed (repeated from Table 2) compared with root- 
mean-square errors of RAOB-derived standard geo- 
strophlc and modified Ekman wind data 82 

5 Mean "actual" values of u, v and wind speed (re- 
peated from Table 2) compared with mean arithmetic 
(bias) errors of RAOB-derived standard geostrophlc 

and modified Ekman wind data 83 

6 Correlation coefficients comparing "actual" wind 

fields with RAOB-derived standard geostrophlc and 
modified Ekman winds 84 

7 Same as Table 2, except for satellite-derived mod- 
ified Ekman winds ^8 

8 Root-mean-square errors and mean arithmetic (bias) 

errors of satellite-derived modified Ekman winds 
(m s~l) • • 

9 Same as Table 6, except comparing satellite- 

derived modified Ekman winds with "actual" and 
RAOB-derived modified Ekman wind fields 96 

10 Weight function constants selected to achieve ap- 
proximately 50% response at the indicated wave- 
lengths using the Barnes (1973) objective analysis 
scheme 


vii 


LIST OF FIGURES 


Figure Page 

1 Schematic examples of how additional downward mo- 
mentum transport, resulting from geostrophlc 
shear In the turbulent boundary layer, modifies 
the surface crossing angle from that expected in 

a barotropic atmosphere (after Hoxit, 1974) 17 

2 Vertical depiction of the AFGWC Boundary Layer 

Model (after Hadeen, 1970) 30 

3 TIROS-N sounding sites at 0921 GMT 19 April 1979. 32 

4 Rawinsonde stations in the AVE-SESAME II Experi- 
ment at 1200 GMT 19 April 1979 35 

5 Synoptic conditions at the surface and 850 mb at 

1200 GMT 19 April 1979 37 

6 Sounding sites exhibiting low-level inversions... 38 

7 Outlines denoting the 17 X 20 grid area in which 

objective analyses of input data were performed 
(solid) , and the inner (dashed) region used in 
qualitative and statistical comparisons 41 

8 Resolution profiles tested, with the solid line 

denoting an approximately 50% response at wave- 
lengths of 1000 km 47 

9 Objectively analyzed fields of input RAOB temper- 
atures ("C) 53 

10 Objectively analyzed fields of input TIROS-N tem- 
peratures (“0 55 

11 Patterns of dewpoint depression based on objec- 

tively analyzed fields of input RAOB dewpoint 
temperatures (®C) 59 

12 Patterns of dewpoint depression based on objec- 

tively analyzed fields of input TIROS-N dewpoint 
temperatures (®C) 62 

13 Objectively analyzed fields of input RAOB pres- 
sures (10^ mb) 64 

14 Objectively analyzed fields of input TIROS-N 

pressures (10^ mb) 66 


viii 


LIST OF FIGURES (CONTINUED) 


Figure 

15 

16 

17 

18 

19 

20 


21 

22 

23 

24 


25 

26 


27 


"Actual" wind fields (a at 50, 300 and 1200 
m AGL 

RAOB-derived standard geostrophic wind fields 
(m s”^) at 50, 300 and 1200 m AGL 

RAOB-derived modified Ekman wind fields (m s"^) 
at 50, 300 and 1200 m AGL 

Satellite-derived modified Ekman wind fields 
(m s“^) at 50, 300 and 1200 m AGL 

Actual (solid) , RAOB (dotted) and satellite 
(dashed) wind vectors at 1600 m AGL superimposed. 

^®*^tical profiles of standard deviations (m s”^) 
of the differences between "actual" wind speeds 
and RAOB-derived modified Ekman wind speeds based 
on varying values of h ; corresponding linear 
correlation coefficients are also shown 

Same as Fig. 20, except using satellite-derived 
modified Ekman wind speeds 

RAOB-generated values of Km (m^ s~^) 

Satellite-generated values of Km (m^ s~^) 

Vertical profiles of standard deviations (ra s~^) 
of the differences between "actual" wind speeds 
and RAOB-derived modified Ekman wind speeds based 
on varying values of Km* corresponding linear 
correlation coefficients are also shown 

Same as Fig. 24, except using satellite-derived 
modified Ekman wind speeds 

Vertical profiles of standard deviations of the 
differences between original and perturbed TIROS- 
N thermal data (*C) and their resulting modified 
Ekman wind speeds (m s”^) 

Fields of horizontal moisture divergence (lO”** g 
kg" 3-^) generated from (a) RAOB- and (b) satel- 
ILte— derived modified Ekman winds and moisture 
data at 300 m AGL 


Page 


69 

74 

77 

89 

99 


102 

103 

105 

106 


108 

109 


111 


114 


ix 


LIST OF FIGURES (CONTINUED) 


Figure 


Page 

28 

Radar summaries for 1135, 1435, 1735, 1935, 2035 
and 2235 GMT 19 April 1979 

116 
1 ^ Q 

29 

Fields of moisture advection (10”** g kg”^ s“^) at 
300 m AGL derived from (a) RAOB and (b) satellite 
data 

30 

Fields of horizontal velocity divergence (10“^ 
s-^) at 300 m AGL derived from (a) RAOB and (b) 
satellite data 

i. ly 
121 

31 

Fields of relative vorticity (10“® s”^) at 300 m 
AGL derived from (a) RAOB and (b) satellite data. 

122 

32 

Fields of vertical motion (cm s“^) at the top of 
the PBL based on (a) RAOB and (b) satellite data. 

123 

33 

Vertical profiles of standard deviations (m s“^) 
of the differences between "actual" wind speeds 
and RAOB-derived modified Ekman wind speeds ob- 
jectively analyzed for an approximately 50% re- 
sponse at wavelengths of 900, 1000 and 1100 km; 
corresponding linear correlation coefficients 
are also shown 

1 Aft 

34 

Same as Fig. 33, except for satellite-derived 
modified Ekman wind speeds 

i. HO 

149 

152 

35 

Data points used in scale analysis 

36 

Stress profiles based on RAOB- (solid) and satel- 
lite-derived (dashed) modified Ekman data, as 
well as "actual" (dotted) winds 

155 


X 


1 . INTRODUCTION 


a. Baekgpound 

Radiosonde obsenratlons (RAOBs) have traditionally been 
the primary source of meteorological data above the earth* s 
surface. Increasingly, however, weather satellite systems are 
becoming a major contributor of valuable sounding data that 
were previously unattainable. As noted by Scoggins ^ al . 
( 1981 ), this new data source can greatly enhance our knowledge 
of atmospheric structure in that; 1) satellite soundings can 
be made on a global scale with greater time and space resolu- 
tion than is possible with the current RAOB network; 2) all 
measurements are made by the same instrument, eliminating 
errors resulting from the variability between radiosondes; anH 
3 } a satellite meakires the entire vertical extent of the 
sounding almost instantaneously, thereby eliminating errors due 
to downstream drift of radiosonde balloons. 

An inherent difference between a satellite sounding system 

and a conventional radiosonde is that a satellite measurement 

describes mean temperature and humidity within a volume of the 

atmosphere, whereas a radiosonde measures these parameters on a 

point- bv- point baala. And, as noted by Shen et al. (1975), 

"although the radiosonde is capable of observing 
with much higher vertical resolution than is the 
remote satellite sounder (meters compared to kilo- 
meters) , the satellite sounding system usually can 
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achieve much higher spatial resolution and more 
meaningful horizontal gradients than those possible 
with a practical radiosonde network*. 

Another difference between the two systems is that the presence 

of clouds adversely affects satellite soundings, but has little 

effect on radiosonde measurements. 

One way in which our knowledge of atmospheric motion is 
being significantly improved is by obtaining wind fields from 
^ata. Current satellite sensors do not directly 
measure winds; however, these systems do provide sufficient 
information from which wind profiles can be computed. Two 
popular techniques for ascertaining wind fields are cloud 
tracking and use of thermal wind relationships. 

Cloud-derived winds are based on the aissumptlon that cer- 
tain clouds, especially smaill cumulus and cirrus, move at or 
near the environmental wind velocity in which they are embed- 
ded. This assumption has been verified through la situ air- 
craft studies by Hasler at al. (1977) and Hasler at al. 
(1979). Several earlier studies found favorable relationships 
between cloud-derived winds and the observed flow, notably 
Hubert and Whitney (1971), Fujlta at at. (1975), and Suchman 
and Martin (1976). More recent studies have employed satellite 
cloud motion techniques to successfully compute divergence 
(Peslen, 1980), kinematic vertical motion (Wilson and Houghton, 
1979) and low-level moisture convergence (Negri and Vender 
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Haar, 1980). Although cloud tracking has been shown to be a 
valuable method for describing wind flow, the procedure has 
several limitations.' First, the time required for cloud track- 
ing at several levels nay be extensive, even when interactive 
computer methods are used (Johnson and Suchman, 1980). Second, 
trackable clouds may not be present at the areas and levels of 
Interest, or may be obscured by clouds at higher levels. 

Application of thermal wind relations to satellite-derived 
temperature data provides an alternate or supplemental source 
of Information about atmospheric flow. This procedure calls 
for the addition of geopotential thickness, derived from satel- 
lite data, to known height values at a "tie-on" level (often 
the surface), thereby yielding the geopotential field at the 
level of Interest. Fr<Ms this field, geostrophlc winds may be 
computed. Favorable results from the thermal wind technique 
have been obtained by exami nlng cross-sections through baro- 
cllnlc zones which were based on Nimbus 5— derived temperature 
daU (Smith and Woolf, 1974; Shen Ai., 1975; Smith 
1975; Arnold atgl., 1976; Horn it ^1. • 1976). Peterson and 
Horn (1977) used Nimbus 6 thermal data to obtain geopotential 
heights and geostrophlc winds; their results showed good agree- 
ment with bracketing data from the National Meteorological 
Center (NMC) and had good continuity during the three day 
period of study. Moyer it al. (1978) determined that Nimbus 
6— derived winds were sufficiently accurate to describe the 
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major features of synoptlc-soale systems apparent on constant- 
pressure charts obtained from RAOB data. 

Several more elaborate techniques for obtaining 
thermally-derived winds have been tested recently using TIROS-M 
soundings. Manouso and Endlleh (1980) employed TIROS-N data In 
OMiblnatlon with GOES (Geostationary Operational Environmental 
Satellite) cloud— motion measurements to test several procedures 
for determining satellite-derived wind profiles. Their tech- 
niques Included (a) using the thermal wind to build upward or 
downward from the cloud-notion level, (b) Introducing the bal- 
ance relationship to (a), and (c) using an eigenvector 
approach. The eigenvector procedure, based on the use of an 
Initial gradient wind profile that Is built upward from the 
given surface winds, provided the most representative results. 

The studies described above were concerned primarily with 
deriving wind fields In the free atmosphere — that part of 
the atmosphere beginning above the ■geostrophlc (or gradient) 
wind level*. To the author's knowledge, the only study thus 
far to apply satellite thermal data to wind computations in the 
planetary boundary layer (PEL) was performed by Carle and Scog- 
gins (1981). Sounding data from Nimbus 6 were used to derive 
geostrophlc wind fields at the mandatory pressure levels (850 
mb, 700 mb, etc.); the logarithmic wind law (Hess, 1959) was 
then applied to "extrapolate" through the PBL (selected as < 
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150 m) to compute surface wind speed and direction. Surface 
winds were not calculated In the western United States due to 
the method's Inability to handle large variabilities In terrain 
height. However » In areas where surface winds were computed, 
the logarithmic wind law produced results which agreed well 
with observed values In most regions. 

b. Statement of the pt»ohl em 

Further study Into the feasibility of using satellite- 
derived thermal data to obtain wind fields In the planetary 
boundary layer seems Justified. With satellite technology pro- 
viding an ever-improving product (the moat recent la the 
geosynchronous VAS* systaa) , It Is reasonable to assume that 
research aimed towards blending satellite data directly Into 
operational boundary layer models will soon begin. Such 
research, If successful, offers several potential benefits. 
Satellite— derived PBL winds could yield an Improved method of 
deter m i n i n g low-level horizontal moisture convergence, already 
shown to be a precursor of many severe storm outbreaks (Hudson, 
1971; Hegrl and Fonder Haar, 1980). Procedures dealing with 
other kinematic parameters (e.g., low-level divergence and vor- 
tlclty) would also benefit from the Improved data source. The 
addition of satellite-derived data Into stability indices which 

• Visible Infrared Spin Scan Radiometer (2ISSH) Atmospheric 
founder 
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incorporate low-level wind speed and direction (e.g., the Air 
Force Global Weather Central (AFGWC) Severe WEAther Threat 
(SWEAT) Index) would help to better define potential areas of 
severe convective weather. 

The research documented in the following pages expands 
upon the initial effort of Carle and Scoggins (1981) by examin- 
ing the feasibility of using satellite-derived temperature data 
to compute winds at several levels within the planetary boun- 
dary layer. The study is diagnostic in nature; prognostic and 
operational applications are subjects for future research. The 
main goal of the investigation is to obtain meaningful horizon- 
tal wind fields and vertical wind profiles in the lowest 1600 m 
of the atmosphere. A secondary goal is to combine these winds 
with humidity data to define areas of horizontal moisture con- 
vergence, hopefully leading to the identification of potential 
regions of Intense convection. 

TIROS-N sounding data collected during the 1979 Atmos- 
pheric Variability Experiment-Severe Environmental Storms and 
Mesoscale Experiments (AVE-SESAME *79) are incorporated into 
the thermal wind relation and modified Ekman equations con- 
tained in the AFGWC Boundary Layer Model (AFGWC-BLM). 
Described in detail by Hadeen and Friend (1972), the model has 
been operational since March 1969 and is based on work by Ger- 
rity (1967). In particular, data from AVE-SESAME II (19 April 
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1979) are employed. RAOB soundings for this period are used as 
a basis for comparing the computed satellite-derived wind 
values.. 

Before outlining In greater detail the methodology of this 
Investigation ( Chapter 5 ) , It Is prudent to review the 
phenomenology of the problem. An examination of the charac- 
teristics and classic formulations of the PBL (Chapter 2) will 
be followed by pertinent details of the AFGWC Boundary Layer 
Model (Chapter 3). In Chapter the TIROS-N and RAOB sounding 
data will be described. Finally, results of the study will be 
presented In Chapter 6. 


2. REVIEW OF THE PEL THEORY 


a. £fiL structure 

The planetary boundary layer occupies the atmospheric 
domain nearest the earth, wherein wind shear and/or thermal 
convection give rise to small-scale turbulence (Deardorff, 
1972). As noted by Sutton (1977)> the PBL may be subdivided 
Into two layers: A shallow surface (constant flux) layer , 
lying Immediately adjacent to the earth, extends to a height of 
'^'50 m and Is characterized by eddy stresses which are nearly 
constant with height and which are an order of magnitude 
greater than the horizontal pressure gradient and Coriolis 
forces. The vertical wind speed profile, with wind direction 
assumed constant with height. Is determined primarily by the 
vertical temperature gradient and the nature of the underlying 
terrain. 

The Bifwan (transition) layer extends from the surface 
layer to the lower boundary of the free atmosphere (where the 
atmosphere Is treated as an Ideal fluid In geostrophlc equili- 
brium), a height of 1-2 km. In this layer of transition from 
the disturbed flow near the surface to the (supposedly) smooth 
flow of the free atmosphere, the wind profile is determined by 
a balance between the pressure gradient force, Coriolis force, 
and residual frictional effects of the earth's surface. Here, 
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turbulent eddy stress is the same order of magnitude as the 
pressure gradient and Coriolis forces. Above the PBL, turbulent 
transport is caused by clear air turbulence (associated with 
internal gravity waves), towering cumulus clouds, and terrain 
effects on a scale large enough to cause upward propagation of 
energy through the PBL (Deardorff, 1972). 


b. Classic surface layer formulations 


Parameterization of turbulent transport has proven histor- 
ically to be the main stumbling block for PBL modelers. The 
eddy viscosity method, or "K- theory", has long been used to 
simulate boundary layer flow (Estoque, 1973; Sutton, 1977; 
Krishna, 1981). Analogous to molecular turbulence theory, eddy 
stresses are assumed to be proportional to the product of the 
shear of a scalar quantity (such as heat or a component of 
momentum) and an eddy exchange coefficient (eddy conductivity 
or eddy viscosity, respectively). In terms of momentum 
transfer and the mean vertical wind shear, for example, 


and 


— mm ~ - -- OIL 

* p UW ^ Kfg^ pr — — 

oz 

‘Tyt = -piTW = K^p 1 ^. 


( 2 - 1 ) 


where (= -ptx*’ti^’) is the appropriate Reynolds* stress, p 
denotes density and is the eddy viscosity. Since the eddy 
stresses are expressed directly in terms of the gradients of 
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mean motion, K-theory is referred to as a "first-order closure" 
scheme. Prandtl (193^) extended the exchange coefficient 
hypothesis by Introducing the nHy<ng length, s imi lar tp the 
molecular mean free path: 




(S-2) 




Where the mixing length { represents the mean distance that a 
turbulent eddy with excess momentum u’(orv') will travel before 
blending In with the environment. In a sm*face layer exhibit- 
ing neutral static stability, the mixing length is zero at the 
surface and Increases linearly with height, such that I = kz 
where k is the (dimensionless) von Karman constant ('^'0.4}. Con- 
sidering the subscripts in (2-1) to be "understood", substitu- 
tion from (2-2) yields the genereJ. expression (see for example. 


Hess, 1959) 


T = (kzfp 

dz 


du _ u. 
dz kz 


( 2 - 3 ) 


where u» = |— | is termed the 

viscosity can thus be rewritten 


The eddy 


4r = • 


( 2 - 4 ) 


Integration of (2-3) produces the well known 




profile for a neutral surface layer 


It* 9 

u(z) = — In-^ , 
k z. 


( 2 - 5 ) 
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where is the roughness parameter (denoting the level at 
which the mean wind speed is presumed to vanish) . Under non- 
neutral conditions, (2-3) may be modified by employing the 
similarity theory of Monin and Obukhov (Haltlner and Williams, 
1 980 ) , such that 


du 

dz 


= ^ ). 


( 2 - 6 ) 


where is an empirical function expressing the adjustment of 
the velocity profile due to conditions of non-neutrsd stability 
( = 1 for neutral stability) , and L is the Monin-Obukhov 
length 

L = where ^ ^ . (2-7) 

QK * •ti 

which indicates the height at which the magnitude of mechanical 
turbulence (due to Reynolds' stresses) equals the magnitude of 
thermally Induced turbulence. Here, T denotes temperature, g 
indicates acceleration due to gravity, and d represents poten- 
tial temperature. Integration of (2-6) in the manner shown by 
Priestly (1959) yields 


u(s) = ^ In 


u» az 


( 2 - 8 ) 


k L ' 

where a = 0.6 is an empirical expansion coefficient. This may 
be termed the lflg-JlngriAn9ar wind profile . From (2-3) and (2- 


6) it can be shown that the eddy viscosity for the non- neutral 
surface layer is 

9>m(^) 


( 2 - 9 ) 


12 


while the mixing length becomes 

kz 


I = 




( 2 - 10 ) 


Similarity theory and numerous empirical studies indicate 
that there are two principal turbulence regimes, viz., free and 
forced convection. The decision as to which regime is applica- 
ble is often based on the value of the Richardson number, 
which may be expressed (Gerrity, 1967) 

I -2 


s -Z. 


ds 

i9 

dz Cp 

dz 


( 2 - 11 ) 

The magnitude of the wind is denoted by s, while c, is the 


can be 


speolflo heat at constant pressure. The ratio 
expressed in terms of Ri (see for example, Haltlner and Willi- 
ams, 1980), and appropriate expressions for may be obtained 
(Estoque, 1963): 


For_ > -Q.Q3 ( fo rced convection! 


yr* = [fc(z +zJ{l^aR,)f 


ds 

dz 


(2-12a) 


For g -0.03 ( free convection) 


/(■„ = C(z + 



1/2 

d-d 

SL 

■d 


dz 


1/2 


(2-12b) 


An average value of a = -3,0 was deduced from wind and tempera- 
ture data collected for the Great Plains Turbulence Program 
(Lettau and Davidson, 1957). Priestly (1959) has shown the 
constant C to be equal to 0.9. Various authors have studied 
the PEL turbulent flux (e.g.. Dyer and Hicks, 1970; Hess ^ 
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> 1981; Rubensteln, 1981) and wind profile data (e.g., Let- 
tau and Davidson, 1957; Priestly, 1959; Boundary Layer Branch, 
1967; Wlerlnga, 1980) from which can be computed. Whlle/r,n 
Is of the order of 1 m^ s"^ In the Ekman layer (Rao and Hasse- 
brock, 1972), Priestly notes that, for example, near the sur- 
face (viz., 1.5 m above grassland) values around 0.25 m s" 
au*e typical In the presence of winds 4 m s”^ and slight to 
moderate Instability. The greater the Instability, the greater 
will be the value of 

c. jiziaaaifi. Ekman lazan formations 

Extension of K-theory to the remainder of the PEL above 

the surface layer has met with the difficulty of determining an 

appropriate form for the mixing length, upon which Ky^ partly 
depends. The approach classically taken (Hess, 1959; Sutton, 
1977; Haltlner and Williams, 1980) Is to require that the eddy 
viscosity be constant with height and equal to the value at the 
top of the surface layer. 

Characterizing the Ekman layer In Its simplest form — 
that of a steady, hydrostatically stable (neutrally buoyant) 
region — the wind profile Is classically obtained through 
manipulation of the horizontal equations of motion (see, for 
example, Estoque, 1973)* Assuming (a) horizontal mean motion; 
(b) that horizontal mean wind shears are small compared to the 
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vertical mean wind shears; and (c) a balance among the 
Coriolis, pressure gradient and residual frictional forces, the 
momentum equations may be written 


and 


1 ^ 

p dz 


+ fv 


J_ 92. 

p dx 


1 illL _ f,, - 1 i2L 
p dz p dy' 


( 2 - 13 ) 


where t* and Xy are defined by (2-1), and / denotes the 
Coriolis parameter. 


Further assumptions leading to the RVwan equations are 
that (d) Kjn is Independent of height and (e) that the large- 
scale pressure gradient and density do not vary significantly 
with height and may be considered constant. This Implies that 


the shearing stress, r, is a function of height only. Finally, 


It is assumed (f) that the geostrophlc wind Is constant with 

gj^ 

height — i- = 0 , with Its components defined as 
dz 


Uy = - 


J:- ^ 

pf dy 


and 


i, 92 
p dx 


Boundary conditions are chosen such that at 


2=0: u = w = 0, 


and at 


Z Xi H: U =Ug, V ■= Vg, 

where H denotes the top of PBL. Finally, assuming, for simpli- 
city, that the geostrophlc flow is entirely zonal (v, =0), the 
wind profile in the Ekman layer may be written (see, for exam- 
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pie, Holton (1979) 

u = tiy [l — e~*"* cos (mz)] 
and 

V ~Ug [e”"“ sin (mz )], 


( 2 - 14 ) 


where m = 


2AL 




Equations (2-14) combine through the rela- 


tionship p = tan 


-1 


to produce the Ekman spiral, where p 
represents the cross-isobaric angle. This relation describes 
the wind as turning clockwise (veering) with elevation through 
the PBL. 


It must be noted that one or more of the assumptions made 
in deriving the classical Ekman spiral equations (2-14) may not 
be satisfied during individual cases. For example, the eddy 
viscosity is often, in reality, variable with height. As noted 
by Hess (1959) » this should not prevent the wind flow from 
turning clockwise with elevation in the PBL, but may cause it 
to turn at a different rate than that described by the 
"spiral*. Another point is that strong warm and cold air 
advectlon will cause the pressure gradient to change rapidly 
with height instead of remaining constant. Hoxit (1974) exam- 
ined the effects of warm and cold air advectlon on geostrophlc 
wind profiles. His findings suggested that the majority of 
observed vertical wind variability for the lowest 2.5 km is 
related to changes in stability, wind speed and/or the direc- 
tion and magnitude of the thermal wind. Themal stratification 
(stability) determined, to a large extent, the vertical dlstri- 
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distribution of turbulence (momentum transport). When geos- 
trophlc shear existed In the PBL, the vertical momentum tran- 
sport was modified from that expected under barotroplc condi- 
tions; this. In turn, modified the cross-isobarlc angle. As 
noted by Hoxlt, and later documented by Krishna (1981), the 
result Is that W 2 urm (cold) air advectlon decreases (Increases) 
the cross-isobarlc angle of the surface wind and Increases 
(reduces) the veering with height of the Ekman layer winds 
(Fig. 1). 

A third factor concerning the classic Ekman spiral equa- 
tions Is the assumption that the geostrophlc wind does not vary 
with height. This assumption Is valid for b 2 u*otropic condi- 
tions; however, as noted by Arya and Wyngaard (1975), the lower 
atmosphere Is normally highly barocllnlc, with geostrophlc 
shears produced by large horizontal temperature gradients asso- 
ciated with various mesoscale systems. Lastly, If one is 
adhering strictly to the assumptions, (2-14) should be applied 
only to purely zonal geostrophlc wind flow conditions. 

Despite these deficiencies, first-order closure (K- theory) 
formulations, because of their mathematical simplicity, have 
been widely used and have provided acceptable results for many 
applications (Bodln, 1980). 


d. Higher- order olasarg models 



Fig. 1. Schematic examples of how additional downward momentum 
transport, resulting from geostrophic shear in the tur- 
bulent boundary layer, modifies the surface crossing 
angle from that expected in a barotropic atmosphere 
Rafter Hoxit, 1974). Vq denotes the surface wind vector; 
Vq(B), the surface wind vector for barotropic conditions; 
^g, the geostrophic wind vector; the surface geo- 

strophic wind vector; Uq, the angle beti^^een the surface 
wind vector and surface isobars; and g, the angle mea- 
sured clockwise from the direction of to the direc- 
tion of the mean thermal wind in the lowest 100 mb. 
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In recent years, various turbulent energy (higher-order 
closure) schemes have been under development. While the 
first-order schemes use the turbulent transfer coefficient K 
for momentum flux), the higher-order closure models use 
transfer coefficients in Integrated flux form (Businger, 1973). 
The result is a considerable complication of the closure prob- 
lem, leading to many additional equations for a complete model 
(see, for example, Mellor and Yamada, 1974; Wyngaard et al . . 
1974; Arya and Wyngaard, 1975; H2atiner and Williams, 1980). 
Such models require fairly fine-resolution input data (e.g., 
the Wangara Experiment data, described by Hess ^. , 1981). 
The satellite data available for this study, however, have a 
much coarser resolution, with only one time period available 
for use. Although the higher-order schemes provide an improve- 
ment upon the classic K-theory formulations as regards detail 
of analysis, such procedures are beyond the scope of the 
current investigation, which focuses on a simpler parameteriza- 
tion of the mean properties of the PEL. 


3. THE AFGWC BOUNDARY LAYER MODEL 


The AFGWC Boundary Layer Model generates prognostic tem- 
perature values which are inserted into diagnostic wind equa- 
tions to yield three-dimensioncLL wind fields in the lowest 1600 
m of the atmosphere (Hadeen and Friend, 1972). Based on work 
by Gerrity (1967)» the horizontal wind flow is deduced using a 
modified version of the Ekman spiral equations described in 
Chapter 2. A goal of the current research is to obtain mean- 
ingful "satellite-derived* horizontal wind fields by inserting 
observed TIROS-N thermal data into Gerrity *s modified Ekman 
equations. Pertinent equations from the model are presented in 
this chapter, with detailed derivations provided in Appendix A. 

a. Fundangntal aspects st Itis. model 

A terrain-following coordinate system derived for a spher- 
ical, rotating earth (Gerrity, 1967) is employed in the AFGWC- 
BLM. The model coordinates are defined 
X = a n{fi) cos <p cos 

y ~ a n{<p) cos ^ sin X, 
s = r, - o - E(x,y). 

where (r,,;»,X) are the spherical coordinates; *a” is the mean 
radius of the earth; £(x,y) is the elevation, in meters, of the 
terrain above mean sea level; and n(t») is the polar stereo- 
graphic map scale factor. 


20 


The model defines the PBL as the region between z = 0 and 
z=H, where H s 1600 m (Hadeen, 1970). The surface layer 
occupies 0 ^ z ^ /i where h = 50 m, while the transition layer 
occupies h ^ z ^ H. 


Based on the aboye-descrlbed coordinate system, the hor- 
izontal equations of motion may be written: 


and 


du 

eU 


dv 

dt 



^ + no ^ 

. JL 

tr 

p 

dx dx 

dz 


__1 

9 

dp dE 

dz 

y dv 
dz; 


(3-1) 


(3-2) 


Here, the following assumptions have been made: (1) all terms 
Inyolving yertical yelocity which are not p«u*t of the indiyl- 
dual deriyatiye are neglected and (2) all terms dealing with 
the yarlablllty of the map scale factor, n(^p), or the conyer— 
gence of meridians are neglected. The map scale factor is 
taken to equal unity, its value at the standard latitude for a 
polar stereographio projection. 


b. Surface layer eouations 

A desirable aspect of the AFGWC-BLM is that it treats the 
PBL In terms of the mean flow, rather than requiring use of 
flux-form data. This circumvents difficulties in defining the 
Monin-Obukhov length, L in (2-7), when flux data au?e not avail- 
able. With the assumption that eddy fluxes in the surface 
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layer are independent of height (noted in Chapter 2, Section 
a), we may write the following relationships (Gerrity, 1967): 


dz 


and 


dT 

dz 




(3-3) 


(3-4) 


where K is the eddy exchange coefficient, s is the horizontal 
wind speed, and is a constant with units of °k. 


Recall (frcmi Chapter 2, Section b) that two types of tur- 
bulence regimes exist: free and forced convection. When 
Ri S -0.03, free convection (buoyancy force) is assumed to dom- 
inate, and the eddy exchange coefficient is given by (Priestly, 
1959) 


K = Xr* 

1. 

IL+jl 


1 ^ 

dz Cp 


When > -0.03, however, forced convection (inertial force) is 
assiimed to prevail and the eddy exchange coefficient is given 
by 

/r = [fcz(l-/J/?J]2|£.. (3-6) 

Gerrity defines the dimensionless quantities X and p to equal 
1.2 and 2.0, respectively. (From this point on, fi no longer 
denotes the cross- isobarlc angle.) Appropriate surface layer 
equations for both convection regimes are given below (see 
Appendix A for derivations): 


Foroed 


regime 


Equations defining wind speed and at z - h during 
forced convection conditions are, respectively: 


JC^k — fc/iu»[l,0 + 


T? u» 


£CS.S convection j:sgjj|£ 


Equations for wind speed and K„ at z = h during free con- 
vection conditions are, respectively: 


c / ..V Zo + l 

5 = s(z=h) = — In — 

A 


Q Q 2 

jL - (z 


A",* = 0.49 X*/i* ^ |u.T>. 


(3-10] 


Expressions for variables common to both the free and 
forced equations are: 

«.*» ( L, , .7 

if* = — — 7-77 r 1.0— ,1.0+ — + -^ I 

2fig k{h-Zi) ^u?[fi-Zi c, J 

( 3 - 11 ; 

u, = C (0.07825 - 0.00825 log R,). (3-12; 

n 

7?o ss — — (Rossby Number). (3-13] 

/*0 

where C represents the magnitude of the surface geostrophlc 
wind, and Tk denotes the temperature at z - h (obtained from 
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sounding data) . 7\ is the teioperature at z - Zi, the height of 

the Instrument shelter; for the purposes of this study, Ti s 
Tgfe and z< = 1 m. 

The angle of deviation, if , between the wind in the sur- 
face layer, S, and the surface geostrophic wind, G, is an empir- 
ical relationship given by Gerrity as 

i> = 0.62S(log /?,)*- 12.750(log /?„) + 80.625, (3-14) 

where if, in degrees, varies between 15.6° for Ro = 10*° and 
32.5° for = 10®. Expressions for the components of S, viz., 
U and /, may be obtained from relationships among S, f and 
C (=[Ci + 

U = (C,iT, - GyKz) id + C„2)-i 
and 

V = {GyKi + CgKz) {Gx + Cy)~*. (3-16) 

where 

JCi = S (C/ + G^Y'^ cos if (3-17) 

and 

Kz^ S (C* + GyY'^ sin if. (3-18) 

c. Computational scheme for the surface layer 

A determination of the appropriate stability regime within 
the surface layer must be made at each computational point; the 
choices for the model are (1) strongly stable, (2) forced con- 
vection, and (3) free convection. As noted previously, /?< is 
often used to indicate the extent of turbulence present. 
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However, the classical definition for Ri (2-11) cannot be used 

ds 

in this research because the vertical wind shear ( — in 2-11 

dz 

and 3-6) is not Inmedlately discernible fr<»i the thermail data. 
Gerrlty has derived an alternate expression for Ri in terms of 
variables already presented, viz., 

Ri = i 1.0- 1.0 + ^ hh (3-19) 

/J I [ ti/ j 

Using this equation, determination of the appropriate convec- 
tion regime is made as follows. 

The first decision to be made is whether conditions are 
strongly stable or forced. The limiting factor for forced con- 
vection on the "stable side" is that must be a positive 
term. Recalling (3-8), -d* is the only variable which may be 

negative. Therefore, can become negative when 

[1.0 + ^ < 0. (3-20) 

lU 

Substituting zero for this term in (3-19) yields 



Thus, when Ri > 0.5. the surface layer may be considered strongly 
stable, and the equations for .S and are 

S= 0.176 C (3-21) 

and 

Km = 1.0 m2 s-‘. (3-22) 

The fraction 0.176 was estimated by Gerrity from the classic 
Ekman spiral in which Ks 1.0 m^ s”\ /= 10“^ s“^ and the wind 
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Is assumed to vanish at z -0, Equation (3-22) Is the classi- 


cal minimum value of 


The determination of whether the surface layer Is In a 
free or forced convection regime must consider the Internal 
constraint In the expression for -d* (3-11). The value for d* 


will be complex If: 

2 

< 0. (3-23) 

Such a situation will occur when the "surface** temperature, Ti, 


1.0 + -JfiL + JL 


h—z^ 


In(^) 

Zi 


Is several degrees warmer than the temperature at s s 50 m, 7\ 


— an unquestionably free convection condition. In other 


words, the equation for d« may be considered valid for forced 
convection conditions as long as (after manipulation of ( 3 - 

23 )): 

Cp J 

Further, as noted previously, turbulence is assumed to be of 
the forced convection type when 




-0.03 < S 0.5. 


(3-25) 


Thus, for this study, the tests for determining the 
appropriate turbulence regime in the surface layer are: 

1) For Jii>0.5, strongly stable conditions are 
assumed, and ( 3 - 21 ) and ( 3 - 22 ) are used to compute 5 and/T„,. 

2) If (3-24) ansi (3-25) are satisfied, forced con- 
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veotlon conditions are assumed to exist, and (3«7) and (3-8) 
are used to compute S' and If the resultant value for Is 
found to be negative, both parameters are recomputed using the 
free convection equations. 

3) When Ri S -0.03, free convection conditions are 
assumed, and (3-9) and (3-10) are used to compute S and A'„i. 

This computational scheme differs In several respects to 
that outlined by Gerrity, most notably In the equalization in 
(3-24). 


d. Icanampn lazac equations 


Recall assumptions (1) and (2) presented in Section a. As 
noted by Gerrity, (3-1) and (3-2) may be further simplified for 
the transition layer by assuming: (3) the Individual deriva- 
tives may be neglected (see Appendix D), (4) Independent 
of height, and (5) the geostrophlc wind components, defined as 


li, = - 


1 

Pf 


dy 


^ dE 


( 3 - 26 ) 


and 


^9 


1 

Pf 


iE. 

dx 


P9 


M. 

dx 


( 3 - 27 ) 


are the following linear functions of height: 


= tt^ + 6 {H—z ) 

( 3 - 28 ) 

Vg = + c(H-z). 

( 3 - 29 ) 
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Substituting (3-26) and (3-27) into (3-1) and (3-2), 

respectively, the horizontal equations of motion may be written 
f ^ 


(recall m = 


-g = -2 m* (v-vg) 


= 2 m* {u-Ug). 


Choosing boundary conditions such that 


V -*Va 


as z -» oo (s /7 for the PBL) 


u -» £/ 
V -* V 


as z h. 


the solutions to (3-30) and (3-31) yield the modified Ekman 


equations I 


u =-Ug + cos[m(z-/i)] + (V-v^) sin[m(z -/i)]{ 


(3-32) 


w = Wg cos[m(z-/i)] - (U-u^) sin[m(z-/i)JJ. 

(3-33) 

Constants 6 and c. in (3-28) and (3-29) are obtained through 
manipulation of the Hypsometric Equation and Equation of State; 


the resulting expressions for Ug and v,, are: 

ti- = 14#+ (H-h) — 1= 

^ ^ H-h Ta * / ^ '' dy t* 


(3-34) 


«. = vf . . 7-. i (K-.) xU 


(3-35) 


where 
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and 


/ 9y’ 


a at' 

V = . 

/ dx 


(3-36) 


(3-37) 


The geostrophlc wind components at the top of the transition 
layer, and are computed using (3-26) and (3-27) • Values 
for tig 'and Wg are obtained using (3-34) and (3-35) for z = /i, 
7a and are temperatures at z = /i and z = H, respectively, 
while = is evaluated over the entire layer, h i z These 
modified geostrophlc wind equations, (3-34) and (3-35), are a 
variation of expressions for the thermal wind (see, for exam- 
ple, Arya and Wyngaard, 1975), and serve as input to the modi- 
fied Elonan equations (3-32) and (3-33). 


e. Computational scheme XSC IhS. transition layer 

Since the surface layer is tcdcen to be the lowest 50 m of 
the PBL, wind flow within that narrow region has not been cal- 
culated in this study. However, data at its boundau*ies (sur- 
face and 50 m) are used to compute the wind components and the 
eddy viscosity coefficients at the upper boundary (lowest level 
of the Ekmaa layer) . The surface parameters are then input 
into the modified Ekman equations to obtain horizontal wind 
fields at 50, 150, 300, 600, 900, 1200, and 1600 m above ground 
level (AGL). The fixed geometrical thickness between levels 
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Insures that a detailed description of the boundary layer winds 
can be obtained over variable terrain (Fig. 2) . 



ajnssajd 
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Fig. 2. Vertical depiction of the AFGWC Boundary Layer Model 
(after Hadeen, 1970). 



4. DATA 


a. Satellite 

TIROS-N sounding data, processed by the National Environ- 
mental Satellite Service (University of Wisconsin) , were made 
available through the auspices of NASA*s Marshall Space Flight 
Center, Huntsville, AL. The data retrieval and processing 
techniques are described In Smith at al. (1978). The slxty- 
two satellite soundings (Fig. 3), obtained at 0921 GMT 19 April 
1979» Included temperature, dewpoint, and geopotentlal height 
data at the surface and at ten pressure levels (1000, 850, 700, 
500, 400, 300, 250, 200, 150, and 100 mb). Smith at (1981) 
have shown that TIROS-N soundings have a nominal global hor- 
izontal resolution and spacing of 250 km, with a higher hor- 
izontal resolution of 50 km being possible for limited areas 
when man-machine Interactive processing methods (e.g., McIDAS) 
are employed. The soundings used In this study have an average 
spacing of approximately 150 km. 

Several studies have evaluated the representativeness of 
TIROS-N temperature data through direct comparisons with 
corresponding RAOB soundings (Phillips ^ ai. , 1979; Smith sL 
1979; Schlatter, 1980 and 1981; Scoggins et al.^ 1981). 
Schlatter (1980 and 1981) found that the average TIROS-N sound- 
ing was too warm neeu? the surface, too cool In the middle tro- 
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posphere, and too warm at the tropopause. Root-mean-square 
(RMS) differences were < 2°C in the 1000-850 mb layer. The 
study also indicated that horizontal temperature gradients 
inferred from TIROS-H soundings were consistently weaker than 
analyzed gradients. A comprehensive study by Scoggins ^ al . 
(1981) compared TIROS-N data with radiosonde soundings during 
AVE-SESAME I (10-11 April 1979). Results showed a mean differ- 
ence of -0.5°C over the entire sounding, Indicating that the 
TIROS-M-derived temperatures contained a negative bias relative 
to radiosonde-derived temperatures. Further, a mean RMS 
difference of 1.8®C was found. 

Cloud contamination is a factor in producing discrepancies - 
between satellite and radiosonde data. Schlatter (198I) noted 
that discrepancies are enhanced in cloudy regions, such that 
cloudy retrievals are more than 1°C too wsirm in the 1000-850 mb 
layer. In the lower troposphere, RMS differences for partly 
cloudy soundings eure laurger than those for clear soundings. 
Hillger and Yonder Haar (1977)» Phillips ^ al. (1979) and 
Smith (1979) also have documented the effect of clouds 

on satellite temperature retrievals. 

The above-mentioned studies provide general guidelines 
about the nature of the input satellite data. However, com- 
parison statistics are not a true measure of the satellite's 
accuracy (Yates, 1974; Bruce ^ aL- , 1977; Phillips sL al- , 
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1979; Smith st al.«» 1979). As noted earlier, discrepancies 
between the two sensors 80*6 due in part to the fact that a 
radiosonde gives a point measurement while the satellite- 
inferred value represents avereige conditions over a volume 
which is several tens of kilometers wide and several kilometers 
in depth. Also, there is usually a time difference between a 
given satellite pass amd corresponding radiosonde soundings. 
Hlllger and yonder Haar (1979) more clearly defined the accu- 
racy of satellite thermal data by performing structure function 
analyses of satellite soundings alone. This procedure elim- 
inates those components of the discrepancies attributable to 
sensors not being co-located in space and time. Their study 
found a noise level of 0.5°C for temperatures at 700, 500, and 
300 mb, a value much smaller than the 1-3°C error usually asso- 
ciated with direct satellite-radiosonde comparisons. 

b. Radiosonde 

Rawinsonde data from AVE-SESAME II, viz., 1200 GMT 19 
April 1979» consisted of standeurd upper-air measurements from 
twenty- three National Weather Service (NWS) stations and six- 
teen special sites (Fig. 4). Although denoted as 1200 GMT 
soundings, most releases actually were made at approximately 
1115 GMT. The NWS stations alone yield synoptic-scale resolu- 
tion through an average station spacing of 400 km, while the 
combination of NWS and special site stations results in 
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subsynoptlc-scale resolution with em average station spacing of 
250 km. These data, also provided by the Marshall Space Flight 
Center, were available In two formats: (1) soundings with data 
at 25 mb Intervals and (2) baroswltch contact data (from which 
the 25 mb soundings were Interpolated) at Intervals of approxi- 
mately 10 mb (100 m). Reduction procedures used to process the 
rawlnsonde data are described by Fuelberg (1974). Further 
details about the AVE-SESAME '79 experiments are given by 
liberty ^ (1979) and Hill jst AL. (1979). Particulars 
concerning the AVE-SESAME II period are provided by Williams ^ 
al* (1980a,b). 

c. arapptlg situation 

At 1200 GMT 19 April 1979 a cold front anchored in West 
Texas extended through low pressure In eastern Colorado, and 
northward across western Nebraska and South Dedcota (Fig. 5), A 
weak warm front stretched from South Dakota southeastward to 
Georgia. A few areas of weak-to-moderate convection were scat- 
tered over southern Texas and along the Mississippi River Val- 
ley. Low-level subsidence and radiation Inversions were pre- 
valent at the various reporting sites (Fig. 6). As the day 
progressed, a low-level tongue of warm, moist air flowed north- 
ward from the Gulf of Mexico, auid the northern end of the front 
began to swing across Nebraska Into eastern Iowa. By 1800 GMT 
major thunderstorm activity began developing over much of 



Fig. 6 


. Sounding sites exhibiting low-level inversions. RAOB 
sites are denoted by circles, satellite soundings by 
triangles . 
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Texas, northwestern Missouri and central Iowa, and by 2100 GMT 
a line of Intense activity had formed through central Nebraska, 
central Kansas and western Oklahoma. A detailed synoptic dls> 
cusslon of the entire AVE>SESAME II period (1200 GMT 19 April- 
1200 GMT 20 April) Is presented by Williams ^ (1980b). 


5. METHODOLOGY 


The main thrust of this study is to obtain wind profiles 
in the PBL using satellite thermal data. The mechanics of the 
research are outlined below, idille remarks on statistical and 
qualitative comparisons performed at the various steps are 
presented in Chapter 6. 

a. Generation sL input data 

Soundings from the 62 TIROS-N and 39 rawinsonde (25 mb- 
interval) observations were plotted manually on Skew T, Log p 
diagrams. Temperature, dewpoint and pressure values for each 
site were linearly Interpolated to the seven "above ground" 
model levels from the plotted soundings. As a check on the 
interpolation procedure, constant height maps of the three 
parameters were plotted and analyzed to Insure continuity of 
the fields. 

b. flb.lSStlY? analyses siL input data 

1 ) Grid structure 

The Barnes (1973) objective analysis scheme was used to 
obtain gridded fields of the various input parameters. A 17x20 
base grid, with a mesh of 127 km, encompassed all TIROS-N and 
RAOB sounding sites (see Fig. 7). The grid was oriented to a 
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Outlines denoting the 17 x 20 grid area in «hich objec 
tive analyses of input data were performed (solid line) 
and the inner region (dashed line) used in qualitative 
and statistical comparisons. Terrain elevations are 
given in hundreds of meters. 
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1:10|000,000 polar stereograpbic map projection such that the 
•upper left corner of the grid was anchored at 45°N 107°W, with 
the lefthand side parallel to the 107^W meridian (see Saucier, 
1955). A subset of the total region, comprising the common 
areas of satellite and rawlnsonde soundings (edso seen in Fig. 
7)t was used to evaluate the computed winds and kinematic 
parameters which are discussed in Chapter 6. This smaller area 
contains 111 grldpolnts; a larger statistical data base would 
be desirable, however 111 is the maximum number of "common" 
points feasible using the given data. Thus, statistics dis- 
cussed in this study should be considered for comparison 
(inspection) purposes only, and not as definitive statements 
for all TIHOS-N data. 

2) Traperature and dewpoint 

Separate 17x20 gridded fields of satellite and RACE tem- 
peratures and dewpoints were generated for each model level. 

3) Pressure 

Pressure fields at the eight levels were obtained by first 
reducing the given station pressures to mean sea level (MSL) 
using the standard pressure-height relation: 

Pmat ~ Psta + 0. 1065(z,*a). ( 5 - 1 ) 

Here, pressures at MSL, p,^, and the station elevation, 
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are In mb, while the station heights, are in meters. The 
Barnes scheme was then applied to yield separate 17x20 grids of 
Pmsi both the RAOB cind satellite data. Using terrain 
elevations provided by AFGWC (shown in Fig. 7), "station pres- 
sures" for each gridpoint were computed by reversing the pres- 
sure-height relation such that: 

Ps/e = Pmsi - 0.1065(2gp). (5-2) 
where Pa/c denotes the surface pressure at a given gridpoint 
and Zgp is the terrain elevation in meters at that gridpoint. 
Pressures for the seven "above ground" levels at each gridpoint 
were then computed using the Hypsometric Equation: 


exp 



1 

CM 

1 



29.27 \ 

T* + 273 

f 


(5-3) 


where p, and are pressure (mb) and height (m), respectively, 
at the lower boundary of the model layer under consideration, 
and p 2 and ^2 correspond to the upper boundary. The mean vir- 
tual tonperature, T* was obtained as an avereige of T* at the 
upper and lower limits of the given layer. Virtual tempera- 
ture, 7^, was computed using the approximate relation 

T* « r + (5-4) 


where T* auid the temperature, T, at the level in question are 


”C; the mixing ratio, w, is given in g kg ^and is defined by 


w — 0.822 


( 103 ). 


(5-5) 


The vapor pressure, e (in mb), was obteiined by inserting the 
appropriate dewpoint temperature, TV. into Lowe’s polynomial 
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algorithm (Lowe, 1977): 

e = o, + r/>(oi + To(cl 2 + To (03 + Tj)(a 4 + Tj)(as + oeT*^ ))))), (5-6) 

where, for To in °C, 

ao = 8.107799961. 

01 = 4.436510521xlO~*. 

02 = 1.428945805x10-2, 

09 = 2.850648471x10-*. 

04 = 3.d31240396xl0-®. 

o# = 2.034080948x10-®. and 
oe f 6.136820929x10-”. 

c. Computation siL modified Rinnfln winds 

Equations (3-34) and (3-35) were solved at each gridpoint 
for all above-gi*ound model levels, using temperatures at the 
top {H) and bottom (A) of the PBL generated in Section b-2. 
Equations (3-32) and (3-33) were then solved (also at each 
gridpoint) . Surface geostrophlc wind components, G* and Gy, 
were calculated at a given point using (3-26) and (3-27), where 
we let Ug s Gg and Vg ^ Gy, 

d. J2l standard geostrophlc winds 

Equations (3-26) and (3-27) were also solved at each model 
level to provide standard geostrophlc winds. These fields 
served as a basis for comparisons outlined in Chapter 6. 
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e. Validation of technique 

In the operational AFGVC-BLM, several variables used In 
computing the modified Ekman winds are obtained from various 
prediction algorithms. In the current study, however, all 
variables are obtained solely from observed data. This funda- 
mental difference in procedure requires a csu*eful validation of 
the computational approach. 

A two-step verification procedure was employed. Twenty- 
six of the rawlnsonde sites provided contact ("actual") winds 
at <\,100 m intervals within the PBL (Fig. M). Values were 
linearly Interpolated to the model levels and then objectively 
analyzed as noted earlier to smooth out the mesoscale features. 
If the RAOB-derlved PBL winds (c<wnputed by inserting AVE-SESAME 
*79 data into the BLM equations) compared favorably with these 
"actual" wind fields, the computational approach was judged to 
be suitable. The second step Involved the satellite data. 
RAOB-derlved modified Ekman winds were used as a basis of com- 
parison with their satellite-derived counterparts. A favorable 
agreement established the utility of using the satellite data 
in the BLM scheme. 

f. Response determination 


It would be desirable to take advantage of the relatively 
small horizontal spacing of the satellite (and RAOB) data to 
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determine not only whether satellite thermal data would yield 
reasonable PBL winds, but also whether such winds would be 
representative of the mesoseale. Two factors, however, do not 
permit this Investigation to strive for such a fine resolution. 
First, the technique chosen (modified Ekmem approach) Is based 
on geostrophlc theory and therefore cannot produce mesoseale 
PBL winds which. In reality, are often ageostrophlc. Second, 

It would be desirable to compare satellite-derived PBL winds to 
actual winds as well as to the RAOB-derlved wind data. Since 
the "actual” wind data are from twenty-six RAOB sites with an 
average spacing- slightly less than 400 km (Fig. 4), mesoseale 
evaluations would not be possible. 

The Barnes (1973) procedure allows one to select percen- 
tages of originally resolved amplitudes, as functions of 
wavelength, that are retained by the analysis. This selection 
can be described by means of response curves or resolution pro- 
files. Ten different choices were Investigated for use in the 
current study (Fig. 8). An ''*50% response was obtained for ampli- 
tudes having wavelengths from 300 to 2000 km (right curve) . The 
goal was to select a response that would reduce "noise” gen- 
erated by the random data errors as well as small scale 
features that could not be consistently resolved by the data 
sets, while simultaneously yielding synoptic-scale resolution. 

To determine the response most compatible with the various 
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data sets, the procedure described in Sections b-2 through e 
were repeated for the ten options. Based on statistical 
evaluations between the thermally<-derived winds* and their 
observed counterparts (see Appendix B), the most appropriate 
response appeared to be at a wavelength of 1000 km (solid 

curve in Fig. 8), All results that follow 2 u*e based on input 
data and observed winds objectively analyzed at this resolu- 
tion. This response isolates larger-scale features of the flow 
which should be quasi-geostrophic in nature, except for the 
effects of friction, emd which should be verifiable with the 
available data. It has the added advantage of removing input 
data errors that would appear as short wave phenomenon. 


g. Cflirotatlm SL parameters 


Gridded fields of horizontal moisture divergence. JiD 2 t 
(where negative divergence Implies moisture convergence) were 
ctxBputed using the standard expression: 


- (w^z). (5-7) 

where, again, u;. is the mlxlog ratio (g kg**^) and denotes the 
velocity vector. Expansion of (5-7) yields 




= w 


du ^ §2^ 
dx dy 


. dtu , .. dw 

+ + V — 

ax ay 


= w (^a * ^2 • (%'uj) 

no 2 = (divergence term )2 + (adveetive terrn )2 . 

As noted by Megri and Vonder Haar (1980), horizontal moisture 


(5-8) 
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convergenc« (negatlTe AfDa) can occur when significant horizon- 
. tal convergence [w (^2 pz) <0] and/or advection of moist air 
< 0 ]are present. Moisture convergence fields for the 
seven above-ground model levels were computed in three ways: 
First, a "control" data set was obtained by solving (5-8) using 
RAOB-derived modified Ekman wind components and values for mix- 
ing ratios derived from RAOB dewpoint temperatures. A second 
set of moisture convergence fields was computed using 
corresponding satellite-derived modified Ekman and mixing ratio 
data. As a test of the "usability" of TIROS-N moisture data, a 
third set of horizontal moisture convergence fields was 
obtained for comparison by combining satellite-derived modified 
Ekman wind components with RAOB-derived mixing ratios. 

To more fully examine the nature of the resultant moisture 
convergence values, fields of horizontal velocity divergence 
and BglgtVtrg advection were analyzed using both the RAOB and 
satellite data seta. Values of horizontal divergence {DIV 2 ) 
were calculated using the standeu*d definition, 

( 5 - 9 ) 

while moisture advection was obtained through the last term of 

« 

( 5 - 8 ). 

Computation of relative vorticlty and vertical 

motion (VM) were additional tests of the "accuracy" of the 
satellite-derived winds. While each of the four previous 
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kinematic parameters were calculated at Individual model lev- 
els, vertical velocities were generated at 1600 m AGL using 
stepwise Integration through the PBL. Assuming zero vertical 
motion at the surface and constant density, the Continuity 
Equation yields: 

Wmo.» i: [-Wi'a(2,-*<..)] (5-10) 

<»2 

where summation Is over the eight model levels and DIV 2 denotes 
the mean horizontal velocity divergence In a given layer. DIVz 
at 50 m AGL wa» assumed to be representative of the divergence 
in the surface layer. Since, as noted by O'Brien (1970), essen- 
tially no correction to kinematic vertical velocities is 
required in the lower atmosphere, current values were not 
adjusted. 


6. RESULTS 


Results from the procedure outlined in Chapter 5 are given 
In the following manner: 

First, grldded fields of temperature, dewpoint and pres- 
sure are presented to compare the Input RAOB smd TIROS-N data. 
Next, RAOB-derlved geostrophlc and modified Ekman winds are 
evaluated with "actual” wind fields to determine the validity 
of the computational technique. The goal Is to show that the 
modified Ekman approach Is not only feasible, but that It 
yields better results than those from the much simpler standard 
geostrophlc wind equations. 

With the RAOB-derlved Ekman winds established as a basis 
for comparison, the satellite-derived Ekman winds are then 
examined. A great similarity between the two Ekman wind sets 
will demonstrate the usefulness of satellite— derived thermal 
data In computing synoptic-scale PBL winds. A sensitivity 
analysis then examines the effects of three variables on the 
computed winds, viz., the depth of the surface layer, the char- 
acter of Kn and the computational scheme's handling of random 
error. In the Input satellite thermal data. Finally, fields of 
kinematic parameters derived using satellite data are compared 
to their RAOB-derlved counterparts. 

a. Comparison sL Incuji. 
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Objectively analyzed fields of input temperature based on 
RAOB and satellite data are shown in Figs. 9 and 10, respec- 
tively. From the surface to 600 m AGL, both versions exhibit 
thermal ridges extending from central Texas into Kansas, and 
from Louisiana northward along the Arkansas-Mlsslssippl border. 
Coldest temperatures are located in the northeastern quadrant 
of the area, over southern Illinois. Based on the satellite 
data, the warm tongue located over the southern Mississippi 
River Valley (above 600 m AGL) extends northwestward into Kan- 
sas and Nebraska. With the RAOB data, however, the warm region 
over the Midwest extends froa western Texas. Furthermore, the 
cold area in the northeast attains more of a trough-like pat- 
tern that stretches into eastern Texas. 

The satellite temperatures, shown in Table 1, tend to be 

cooler than their RAOB counterparts. This is observed for 
both the horizontal mean and wa-ri minn values. Greatest aversige 
differences of 3.2 and 3.3°C occur at 600 and 900 m AGL, 
respectively. Standard deviations for the satellite data are 
only slightly smaller than those for the RAOB data. Thus, the 
two sources Indicate similar* horizontal variability even though 
their mean values are slightly different. A portion of the 
differences is assumed to be due to the TIROS-N soundings being 
taken approximately 2 h earlier than the rawlnsonde soundings. 
Nonetheless, the statistical comparisons are in good agreement 
with those from previous evaluations (e.g., Arnold ei. al. , 



Fig. 9. Objectively analyzed f: 

rc). 


of input RAOB temperatures 
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Fig. 9. Continued. 
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Fig. 10. Objectively analyzed fields of input TIROS-1 
tures (®C) . 


tempera- 


55 





Fig. 10. Continued. 
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Table 1. Average and maxlsaum values, as well as standard devia- 
tions (o), of the gridded RAOB and satellite tempera- 
ture, dewpoint and pressure fields. 


RAOB 


Temperature (®C) 

Dewpoint (“O 

Pressure 

(mb) 

Level (m) 

Avg 

Max 

a 

Avg 

Max a 

Avg 

Max 

a 

1600 

11.1 

13.8 

1.4 

3.8 

10.0 6.1 

824 

844 

13.4 

1200 

13.4 

15.8 

1.1 

7.1 

12.8 4.0 

864 

885 

14.1 

900 

15.0 

17.5 

1.3 

7.2 

15.0 8.5 

896 

917 

14.6 

600 

16.0 

19.5 

1.7 

10.0 

16.9 5.4 

928 

950 

15.0 

300 

16.4 

20.8 

2.0 

11.2 

18.8 5.4 

961 

984 

15.4 

150 

16.2 

20.3 

2.0 

12.2 

19.0 4.5 

979 

1000 

15.7 

50 

16.0 

19.4 

2.2 

12.8 

18.3 4.0 

990 

1010 

15.9 

Sfc 

15.7 

19.8 

2.4 

13.0 

18.2 3.7 

996 

1020 

15.9 





TIROS- 

•N 






(°C) 

Dewpoint (“C) 

Pressure 


Level (m) 

Avg 

Max 

a 

Avg 

Max a 

Avg 

Max 

a 

1600 

8.9 

12.2 

1.2 

-1.2 

10.3 4.5 

822 

843 

13.2 

1200 

10.6 

13.7 

1.1 

2.8 

12.3 4.6 

863 

884 

13.8 

900 

11.7 

14.9 

1.2 

5.4 

13.6 4.4 

895 

917 

14.2 

600 

12.8 

16.1 

1.5 

7.7 

15.1 4.1 

927 

950 

14.7 

300 

13.9 

17.3 

1.8 

9.9 

16.4 3.8 

961 

984 

15.2 

150 

14.5 

17.8 

2.0 

10.9 

17.0 3.7 

978 

1000 

15.5 

50 

14.9 

18.2 

2.1 

11.6 

17.3 3.6 

990 

1010 

15.7 

Sfc 

15.1 

18.4 

2.2 

12.0 

17.7 3.6 

996 

1020 

15,7 
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1976; Moyer fit 1978), 

Statistics for the gridded deNpoint temperatures are also 
presented in Table 1. Average values Indicate that TIROS-N 
diagnosed the atmosphere as being drier than when "observed” by 
rawinsondes, particularly at 1200 and 1600 m AGL. Maximum 
values of satellite-derived dewpoints are also consistently 
"drier" (smaller) than their RAOB counterparts, except at 1600 
m AGL a slightly larger value is noted. Based on the 
computed standau*d deviations, less vertical variation is found 
in the satellite data than in the RAOB data. This is indica- 
tive of the volumetric sampling of the satellite's radiometers 
that results in considerably less vertical resolution than is 
obtainable from radiosondes. Also, the satellite standard 
deviations are substantially smaller than rawlnsonde-derlved 
versions in the midlevels, particularly at 900 m AGL idiere 
values are 4,4 and 8.5®C, respectively. 

The relative inability of TIROS-N to capture features of 
the moisture field is clearly revealed in patterns of dewpoint 
depression. In the RAOB data (Fig. 11), the highest moisture 
content is found along a line stretching from Texas to 
Nebraska, This feature is evident at all levels, tilting to 
the east with height. Drier air is indicated along the Missis- 
sippi River Valley, and, above 900 m AGL, in the western third 
of the area. Although differences are evident between the pat- 



Fig. 11. Patterns of dewpoint depression based on objectively 

analyzed fields of input BAOB dewpoint temperatures (**C). 
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Fig. 11. Continued. 
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terns derived from RAOB and satellite data» common features 
Include the southern extent of the moist tongue through the 
midsection of the region, as well as the drier eurea along the 
central Mississippi River Valley (Fig. 12). 

Input fields of pressure from the RAOB and satellite data 
are nearly Identical at all model levels. This can be seen 
both In the statistics (Table 1) and In the horizontal fields 
shown In Figs. 13 end 14. Generally, a ridge of high pressure 
stretches from the Gulf of Mexico to Illinois. A weak trough 
along the eastern Kansas— Oklahmia border Is oriented toward 
southwestern Missouri. Additionally, low pressure extends along 
the western boundary of the region. 

b. Evaluation RAOB-derlved winds 

Fields of horizontal winds obtained fr<xn rawlnsonde ther- 
mal data will now be compared with observed ("actual") pat- 
terns. Based on various statistics that are described later, 
the "best" and "worst" results generally occur at 300 m and 
1200 m, respectively. Thus, gridded vector depictions of the 
various winds will be presented for these levels. In addition, 
fields at 50 m AGL are included to show details at the surface 
layer/transition layer Interface. 

1) "Actual" winds 


"Actual" wind fields at 50, 300, 1200 m AGL are shown in 


Fig. 12. Patterns of dewpoint depression based on objectively 

analyzed fields of input TIROS— N dewpoint temperatures 

rc). 
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Fig. 12. Continued. 
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Fig. 13. Continued. 
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Fig. 14. Objectively analyzed fields of input TIROSH 
(10^ mb). 


pressures 
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Fig. 14. Continued 
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Fig. 15. Meridional flow is pred<MBlnant over the area, with 
the wind dlreotlona veering from south^southeast to south- 
southwest with height. Comparison with Fig. 13 reveals that 
winds In the lowest levels are crossing the Isobars toward 
lower pressure, as Is expected near the surface. The southerly 
flow pattern Is also evident In the horizontal mean values 
given in Table 2. The vertically averaged u-component is -1.9 
m s \ while the average v-component Is much stronger (10.3 n 
s ^ ) . Standard deviations of the v-oomponents su*e greater than 
those for the u-oomponents; vertical aversiges eu*e 3.7 and 2.8 m 
s \ respectively. Vlnd speeds through the PEL average 11.8 m 
, with strongest mean values oocvu*rlng at 600 m and 900 m 
AGL. Standard deviations of the wind speeds at Individual lev- 
els fall between corresponding values for the two components. 

Since "actual” fields are based on only 26 observation 
sites with an average spacing of slightly less than 400 km 
(Fig. 4), the resulting patterns should contain somewhat less 
detail than obtained from thermally-derived winds based on 
either the 39 RAOB soundings (with '^250 km spaclngs) or the 62 
satellite soundings (with '^150 km spaclngs). This result Is 
expected even though each of the three data sets was objec- 
tively analyzed to achieve approximately a 50 % response for 
wavelengths of 1000 km. 


2) RAOB-derlved winds 
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Fig. 15. "Actual” wind fields (m s“l) at 50 . 300 and 1200 m AGL. 
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Fig. 15. Continued. 
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Table 2. Mean values and standard deviations (m s~l) for "actual" 
wind data . Values at each above-ground model level , as 
well as vertical averages, are given for the u- and v- 
conq)onents and wind speed. 



Level (m) 

Mean 

g 

u-component : 

1600 

1.9 

2.9 


1200 

0.8 

3.4 


900 

0.3 

3.0 


600 

-1.2 

2.6 


300 

-4.7 

2.3 


150 

-5.3 

2.6 


50 

-5.1 

2.5 


Average 

1 

• 

VO 

2.8 

v-component ; 

1600 

9.3 

3.8 


1200 

11.4 

5.3 


900 

12.3 

5.4 


600 

12.4 

4.8 


300 

10.2 

3.3 


150 

8.6 

3.1 


50 

7.6 

3.2 


Average 

10.3 

3.7 

Wind Speed: 

1600 

10.4 

3.6 


1200 

12.4 

5.2 


900 

13.2 

5.3 


600 

13.3 

4.5 


300 

12.1 

2.8 


150 

11.1 

2.6 


50 

10.1 

2.7 


Average 

11.8 

3.8 
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Standard geostrophle wind fields, computed using (3-26) 
and (3**27)f are presented as an often-used procedure with which 
to assess the usefulness of the modified Ricwan technique. 
RAOB-derlved geostrophle winds, shown in Fig. 16, appear 
unrealistic at 50 m AGL. Speeds are as great as 33.0 m s~\ 
much stronger than "actual” values (Fig. 15). Differences in 
wind direction are also quite pronounced. The area of light 
and variable winds in the eastern portion of the area is quite 
reasonable since it is located near the center of the hig h 
pressure ridge (Fig. 13). One .should note that this feature 
was not resolved by the "actual” winds due to Inadequate con- 
tact data in the area (Pig. 4). Geostrophle winds at 300 m AGL 
are more compatible with "actual” values, however they still 
are stronger than observed speeds. At 1200 AGL, the RAOB- 
derlved geostrophle winds most closely resemble their "actual” 
counterparts. This improvement with increasing altitude is 
expected since the effects of friction decrease with height. 

Application of the modified Ekman scheme (Chapter 5, Sec- 
tion c) results in the RAOB-derlved winds shown in Fig. 17. 
The Ek ma n winds appear quite comparable to corresponding 
"actual” values (Fig. 15), particularly at 300 m AGL. The 
modified wind patterns appear more reasonable than "actual” 
values in two areas. First, the Ekman fields exhibit a more 
realistic light and variable character in the lower levels near 
the high pressure ridge. Second, along the western border of 



Fig. 16. KAOB-derived standard geostrophic wind fields (m s"l) at 
50, 300 and 1200 m AGL. 
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Fig. 16. Continued. 
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Fig. 17. SAOB-derived mbdlfled Ekman wind fields (m s”^) at 50, 
300 and 1200 m AGL. 
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Fig. 17. Continued. 
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Flg« 17. Continued. 
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the region above 300 m AGL, the Ekman winds depict more coun- 
terclockwise flow, reflecting the presence of a north-south 
oriented trough located Just outside the grid area. Cnnparlng 
Figs. 16 and 17 > It la apparent that the Ekman approach has 
produced mere realistic wind speeds at the lowest levels than 
the geostrophlc technique. Although both geostrophlc and Ekman 
winds reflect the western trough and eastern pressure ridge ^ 
the Rkman vectors flow more toward low pressure than do their 
stronger geostrophlc counterparts. 

Various statistical ccmparlsons between the "actual” and 
two RAOB-derived wind sets are given In Tables 4 through 6. 
First, however, mean values and standard deviations of the 
fields represented by Figs. 16 and 17 are presented in Table 3. 
Horizontal means again reflect fairly light u-components, with 
the geostrophlc values being smaller than RVnian values. Since 
Isobars are predominantly oriented along a north-south direc- 
tion (Fig. 13), the greater Ekman u-components reflect an 
enhanced cross— Isobarlc flow than occurs with observed values 
(Table 2). The Ekman wind speeds 2 ure, correspondingly, smaller 
than their geostrophlc counterparts. Maximum Ekman winds occur 
at 600—900 m AGL (as do the "actual" winds), as opposed to the 
geostrophlc values which continually decrease with altitude. 
Standard deviations for the Ekman data generally are much less 
than those for the geostrophlc winds, partlculeu?ly in the 
lowest levels. The only exceptions are for the u-component; 


Table 3. Same as Table 2 , except for RAOB-derived standard geo- 
strophlc and modified Ekman winds. 
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Table 4. Standard deviations of "actual” u, v and wind speed 

(repeated from Table 2) compared with root-mean-square 
errors of RAOB-derived standard geostrophic 2 ind modified 
Ekman wind data. Units are m s~^. 



Level (m) 

"Actual" 

a 

RAOB— Geo 
BMSE 

RAOB-Ekman 

RMSE 

u-component : 

1600 

2.9 

3.4 

3.4 


1200 

3.4 

3.2 

3.5 


900 

3.0 

2.5 

3.1 


600 

2.6 

2.9 

2.6 


300 

2.3 

5.8 

2.0 


150 

2.6 

6.4 

2.5 


50 

2.5 

6.2 

3.2 


Average 

2.8 

4-3 

2.9 

v-component : 

1600 

3.8 

3.6 

3.6 


1200 

5.3 

3.0 

4.4 


900 

5.4 

2.6 

4.4 


600 

4.8 

3.7 

3.3 


300 

3.3 

8.1 

1.5 


150 

3.1 

10.7 

2.2 


50 

3.2 

12.4 

2.8 


Average 

3.7 

6.3 

3.2 

Wind Speed: 

1600 

3.6 

3.4 

3.3 


1200 

5.2 

3.2 

4.4 


900 

5.3 

2.7 

4.4 


600 

4.5 

3.7 

3.3 


300 

2.8 

7.7 

1.9 


150 

2.6 

10.0 

2.6 


50 

2.7 

11.5 

3.2 


Average 

3.8 

6.0 

3.3 
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Table 5. Mean "actual" values of u» v and wind speed (repeated 
from Table 2) compared with mean arithmetic (bias) er- 
rors of RAOB-derlved standard geostrophlc and modified 
Eknan wind data. A positive bias Indicates that RAOB 
values are greater than "actual" values. Units are m 
s-1. 



"Actual" 

RAOB— Geo 

RAOB-Ekman 

Level (m) 

Mean 

Bias 

Bias 

u-component : 1600 

1.9 

-1.5 

-1.4 

1200 

0.8 

-0.1 

-1.0 

900 

0.3 

0.5 

-1.4 

600 

-1.2 

2.1 

-1.1 

300 

-4.7 

5.5 

1.0 

150 

-5.3 

6.1 

1.2 

50 

-5.1 

5.8 

1.3 

Average 

-1.9 

2.6 

-0.2 

v-component : 1600 

9.3 

0.4 

0.6 

1200 

11.4 

0.2 

-1.0 

900 

12.3 

0.8 

-1.7 

600 

12.4 

2.1 

-1.7 

300 

10.2 

5.8 

-0.4 

150 

8.6 

8.2 

0.0 

50 

7.6 

9.9 

-0.1 

Average 

10.3 

3.9 

—0 .6 

Wind Speed: 1600 

10.4 

0.1 

0.4 

1200 

12.4 

-0.1 

-1.1 

900 

13.2 

0.6 

-1.6 

600 

13.3 

2.0 

-1.6 

300 

12.1 

4.8 

-0.8 

150 

11.1 

6.7 

-0.7 

50 

10.1 

8.3 

-0.8 

Average 

11.8 

3.2 

-0.9 
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Table 6« Correlation coefficients comparing "actual" wind fields 
with R^B-derived standard geostrophic and modified Ek- 
man winds. Linear correlations examine wind speeds, 
while vector values denote total wind comparisons. 


Linear 

Correlation Coefficient* r 

Level (m) 

"Actual" /R-Geo 

"Actual "/R-Ekman 

1600 

0.766 

0.792 

1200 

0.873 

0.764 

900 

0.926 

0.804 

600 

0.937 

0.888 

300 

0.821 

0.931 

150 

0.693 

0.845 

50 

0.696 

0.786 

Vector 

Correlation Coefficient, R 

Level (m) 

"Actual"/R-Geo 

"Actual " /R-Ekman 

1600 

0.642 

0.664 

1200 

0.789 

0.694 

900 

0.859 

0.763 

600 

0.877 

0.850 

300 

0.837 

0.882 

150 

0.763 

0.833 

50 

0.744 

0.801 
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however, one should recall that this component is relatively 
small and variable since the flow Is mainly meridional. Stan- 
dard deviations of the Ekman winds agree better with those of 
the "actual" data (Table 2) than do geostrophic values. Thus, 
the Ekma n procedure better describes both vertical and horizon- 
tal variations in the observed flow than does the geostrophic 
technique. 

Agreements between computed and observed winds can be 
investigated further using Table Here, root-mean— square 
errors of the geostrophic and modified Ekman winds are compared 
with standard deviations of the "actual" data. Examination of 
vertically averaged RMS errors shows that the Ekman winds yield 
smaller values than do their geostrophic count erpeu'ts. The 
same relationship generally holds also for the horizontal wAang 
at the three levels below 900 m AGL. This pattern of Ekman 
winds yielding "better" statistical values than geostrophic 
winds in the lower levels will be reflected in other statistics 
to bo examined as well. It is eissumed to be due primarily to 
handling of frictional effects by the Ekman procedure. The 
standard geostrophic wind equations do not account for fric- 
tion, thereby "degrading* winds in the lower levels. On the 
other hand, the Ekman scheme retains the effects of friction 
equally at all levels (through "degrading", to a slight 
degree, winds in the upper levels. 
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An important point is that Ekman winds have RMS errors 
which are almost always smaller than standard deviations of the 
"actual* wind data (Table 4), Conversely, geostrophic winds in 
the three lowest levels possess RMS errors that are consider- 
ably greater than the observed horizontal variability. The 
Ekmah technique is thus considered the preferred method for 
computing PBL winds. 

In Table 5, systematic errors of the two computational 
methods are examined by comparing mean arithmetic (bias) errors 
of the geostrophic and modified Ekman winds to mean values of 
the "actual" wind data. Comparison of the two rlghthand 
columns shows that a strong positive bias generally results 
from using the geostrophic approach. On the other hand, the 
Ekmw method yields a much smaller negative bias. This can be 
seen particularly in the lower levels where, for example, 
"actual" wind speeds at 50 m AGL are overestimated by 8.3 m s”^ 
using standard geostrophic equations, but are underestimated by 
less than 1.0 m s"^ using the Ekman approach. Although both 
methods yield systematic errors that are smaller than observed 
means (except for u-components in the upper levels), the gtmian 
technique usually produces the least bias. 

Finally, the correlation between the RAOB-derived wind 
sets and the "actual* wind fields is examined (Table 6). Two 
types of correlation coefficients have been employed. The 
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"product moment” formula for the linear correlation coeffi- 
cient, r, Is used to quantify a linear relationship between 
wind speeds of the two data sets. A second pau'ameter, Court's 
total vector correlation coefficient, R, expands upon the 
simpler linear version by including a multiple correlation 
between vectors (Lambeth, 1966). This results in a statistic 
which relates two wind vectors that are separated in time, 
space, or both time and space. The mechanics of both statis- 
tics are presented in Appendix C. 

Examination of the linear correlation coefficients (Table 
6) shows that geostrophlc wind speeds exhibit slightly better 
agreement with their "actual" counterparts in the midlevels 
than do the Ekman data. However, at 1600 m AGL and below 600 m 
AGL, the Ekman wind speeds are more highly correlated to the 
"actual” data, with linear values ranging from 0.786 at 50 m 
AGL to 0.931 3t 300 m AGL. Similar results occur with the vec- 
tor coefficients. Frwn 600-1200 m AGL, the geostrophlc winds 
show slightly better agreement with observed values than the 
Ekman fields (maximum vector value, 0.877). However, at 1600 m 
and below 600 m AGL, the Ekman winds are not only more highly 
correlated to the "actual" winds, but they exhibit the largest 
vector correlation coefficient of all levels (0.882 at 300 m 
AGL). A consideration here is that 300 m AGL is the level at 
which input data are most affected by low-level inversions 
(Fig- 6). Since Inversion conditions more ne 2 U'ly reflect the 
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hydrostatically stable regime upon which the Ekman equations 
are based, data at 300 m AGL would be assumed to yield the most 
representative winds. 

3) Validation 

Qualitative comparisons (Figs. 15, 16 and 17) have shown 
that the modified Ekman technique, coupled with standard RAOB 
data, yields wind fields which agree very favorably to 
corresponding "actual” values. Furthermore, the statistical 
c<Maparlsons (Tables 3 through 6) Indicate that the Ekman 
approach Is not only well founded, but Is preferred over the 
simpler geostrophlc method. Best agreements between Ekman and 
observed winds occur near the surface. These comparisons 
establish the modified Ekman procedure with rawlnsonde input 
as a standard against which to compare satellite-derived Ekman 
winds. 

c. Satellite-derived modified Ekman winds 

Fields of modified Ekman winds obtained from TIROS-N ther- 
mal data are presented in Fig. 18. The satellite-derived winds 
compare quite favorably to the "actual" wind fields in Fig. 15. 
As with the RAOB-derlved Ekman winds, patterns exhibit more 
detail than the observed data by indicating the presence of the 
pressure ridge In the east and. In the upper levels, the trough 
just west of the grid. Mean values and standsurd deviations for 
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the satellite-derived u,v components and wind speeds are given 
in Table 7« Comparison with corresponding "actual" data in 
Table 2 Identifies the satellite winds as having a more wes- 
terly (more negative "u") ccnnponent above 300 m AGL. Vertical 
averages of the v-components show satellite-derived values to 
be '''1.7 m s"^ less than the observed data, resulting in mean 
wind speeds 1,8 m s"^ slower than "actual" winds. The greatest 
mean speeds again occur at 600 m AGL* Vertical averages of the 
standard deviations are quite slollau?, with the maximum differ- 
ence at an individual level being 1.3 m s**^. 

A qualitative comparison of the RAOB- and satellite- 
derived modified Ekman winds (Figs. 17 and 18, respectively) 
indicates a very good correspondence, particularly • in the 
lowest levels. Statistical comparisons between mean values in 
Tables 3 and 7 show satellite Ekman winds have a slightly 
stronger u-component and are ''/0.9 m s“^ weaker than RAOB winds. 
Standard deviations Indicate a somewhat greater degree of vari- 
ability in the satellite data, with differences in RAOB-derived 
values being less than 1,6 m s”^ at any given level. 

Root-mean-square errors between the satellite-derived and 
"actual" wind fields are given in Table 8, Comparisons with the 
two rlghthand columns of Table 4 show that near the surface and 
in the vertical, aversige the satellite errors are between 
corresponding statistics for the RAOB-derived geostrophlc and 


Table 7. Same as Table 
Ekman winds. 

2, except 

for satellite- 

-derived modified 



Level (m) 

Mean 

g 

u-component : 

1600 

0.1 

2.7 


1200 

-0.7 

2.8 


900 

-1.2 

3.0 


600 

-2.7 

2.9 


300 

-3.8 

2.5 


150 

-4.0 

2.1 


50 

-3.7 

2.3 


Average 

-2.3 

2.6 

V“Component : 

1600 

8.6 

5.1 


1200 

9.2 

4.7 


900 

9.6 

4.3 


600 

9.6 

3.6 


300 

8.7 

3.4 


150 

7.7 

3.8 


50 

6.8 

4.1 


Average 

8.6 

4.1 

Wind Speed: 

1600 

9.8 

4.5 


1200 

10.3 

4.3 


900 

10.7 

4.0 


600 

10.9 

3.5 


300 

10.4 

3.0 


150 

9.5 

3.3 


50 

8.7 

3.7 


Average 

10.0 

3.8 


! 





Table 8. Root-mean-square errors and mean arithmetic (bias) er- 
rors of satellite-derived modified Ekman winds (m s"l) . 




SAT-Ekman 

SAT-Ekman 


Level (m) 

BMSE 

Bias 

u-component : 

1600 

4.3 

-1.9 


1200 

4.3 

-1.5 


900 

4.3 

-1.9 


600 

3.7 

■ -1.5 


300 

2.8 

0.9 


150 

3.4 

1.4 


50 

3.9 

1.5 


Average 

3.8 

-0.4 

v-component : 

1600 

5.0 

'-0.6 


1200 

5.8 

-2.2 


900 

5.8 

-2.8 


600 

4.9 

-2.9 


300 

3.0 

-1.6 


150 

2.8 

-1.0 


50 

2.8 

• 

0 

1 


Average 

4.3 

-1.7 

Wind Speed: 

1600 

4.6 

-0.6 


1200 

5.6 

-2.1 


900 

5.6 

-2.5 


600 

4.6 

-2.4 


300 

3.0 

-1.7 


150 

3.2 

-1.6 


50 

3.3 

-1.3 


Average 

4.3 

-1.7 
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modified Ekman winds. Furtheraore, satellite values at Indivi- 
dual levels are no more than 1.6 m s**^ greater than those for 
the RAOB Ekman winds. Smallest satellite-derived RMS errors 
occur In the lower levels, partlcul 2 u?ly at 300 m AGL. The 
satellite-derived Ekman winds have RMS errors that are slightly 
greater than the standard deviations of the "actual” wind data 
(Table 2). 

Mean aurlthmetlc errors of the satellite Ekman winds are 
also presented In Table 8. CcHiparlson with corresponding 
RAOB-derlved errors In Table 5 reveals a slmlleu* pattern of 
negative biases for u-components above 300 m AGL, and v- 
components and winds speeds at nearly all levels. Although 
vertical averages of satellite bias errors are slightly greater 
than those for RAOB Ekman winds, they are significantly smaller 
than those of corresponding RAOB-derlved geostrophlc values. 
Furthermore, these bias errors are significantly smaller than 
the degree of natural variability seen in the "actual" wind 
data (Table 2). 

Correlation coefficients relating the satellite Ekman 
fields to observed winds are presented In the center column of 
Table 9. Values of the linear coefficient show that 
satellite-derived wind speeds are most highly correlated to 
their "actual" counterparts in the lower levels in general, and 
at 300 m ACL In particular (where r equals 0.844). This pat- 


Table 9. Same as Table 6, except comparing satellite-derived mod- 
ified Ekman winds with "actual" and RAOB-derived modi- 
fied Ekman wind fields. 


Linear 

Correlation 

Coefficient, 

r 

Level (m) 

"Actual " / S-Ekman 

R-Ekman/ S-Ekman 

1600 

0.627 



0.842 

1200 

0.632 



0.857 

900 

0.662 



0.870 

600 

0.746 



0.890 

300 

0.844 



0.921 

150 

0.795 



0.941 

50 

0.757 



0.979 

Vector 

Correlation 

Coefficient, 

R 

Level (m) 

”Ac tual ” / S-Ekioan 

R-Ekman/S-Ekman 

1600 

0.334 



0.725 

1200 

0.401 



0.751 

900 

0.498 



0.768 

600 

0.597 



0.784 

300 

0.718 



0.828 

150 

0.735 



0.893 

50 

0.749 



0.957 
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tern is similar to that shown by the RAOB-derived coefficients 
in Table 6. Satellite values are slightly lower than those 
obtained using RAOB Ekman data, yet remain significantly high 
at all levels. 

The added effect of wind direction is evaluated by the 
vector correlation coefficients given in the lower center 
column of Table 9* Note that values below 600 m are quite 
high. In the upper levels, values of R indicate, at first 
glance, a disturbingly low correlation between satellite- 
derived and "actual* winds* Investigation of this discrepancy 
revesiled that the problem lies primarily in comparing two 
"total" wind fields which differ significantly in the "smooth- 
ness" of their flow patterns. As noted previously, "actual" 
fields are based on 26 observations, while RAOB and satellite 
grids are derived fr<mi 39 and 62 soundings, respectively. The 
Increased amount of detail available from the larger RAOB and 
satellite data bases has been shown to be most evident along 
the eastern and western borders of the grid area. Values of R 
which do not consider these border gridpoints reveal a much 
higher correlation between satellite— derived and observed 
winds, raising the lowest value (at 1600 m AGL) from 0.334 to 
over 0.620. 

This problem is not as obvious in the RAOB-derived R - 
values of Table 6 . Superimposing the three types of fields 
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( "actual" » RAOB and satellite) at any given upper level reveals 
that RAOB-derived vectors generally lie between those of 
corresponding satellite-derived and observed winds. Figure 19 
shows, for example, fields at 1600 m AGL. In other words, 
satellite and RAOB vectors exhibit a slmllau* amount of detail 
(Figs. 17 and 18); however, RAOB winds show less directional 
(and speed) departures from observed vectors than do satellite 
winds. 


A more appropriate statistical test of the significance of 
satellite-derived PBL winds, then, is to compare them directly 
with their (more equally detailed) RAOB counterparts. The 
rlghthand column of Table 9 gives the results of such a com- 
parison. Immediately, we see the very high correlation between 
wind si>eeds at all levels, ranging from an r- value of 0.979 at 
50 m AGL to 0.842 at 1600 m AGL. A similar pattern is seen in 
the values of R, which range from 0.957 to 0.725 at 50 m and 
1600 m AGL, respectively. The slight decrease in correlation 
with height probably occurs because the satellite winds retain 
(with increasing altitude) the light and variable pattern along 
the eastern border of the grid, whereas the RAOB winds yield 
stronger and more southerly vectors at the higher levels. This 
may be a result of the denser satellite coverage (recall Figs. 
3 and 4) in the area of the high pressure ridge. 


Synoptic-scale PBL winds computed using satellite-derived 



Fig. 19. "Actual" (solid), RAOB (dotted) and satellite (dashed) 
wind vectors at 1600 m AGL superimposed . 
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temperatures have thus been shown to be very comparable to 
winds obtained in a similar manner using RAOB thermal data. 
Overall, the most favorable comparisons, both qualitative and 
statistical, are at 300 m AGL (the average height of the 
numerous inversions). Satellite winds tend to be slightly 
weaker and exhibit more cross- Isobarlc • flow. Both RAOB and 
satellite modified Ekman winds reproduce the general flow pat- 
tern of the observed fields; given a denser observed data base, 
however, much higher correlation to "actual" winds is antici- 
pated. 

d. Sensitivity analysis 

The sensitivity of three parameters which variously affect 
nearly all terms in the modified Ekman equations is investi- 
gated in this section. Variations in the height of the surface 
layer/ transition layer Interface are examined first. Second, 
the preference for a horizontally variable, rather than a con- 
stant .iTm Is discussed. Finally, the ability of the computational 
scheme to handle random error in satellite-derived temperature 
data is explored. 

Ekman winds are quite sensitive to the appropriate selec- 
tion of the depth of the surface layer because h determines, in 
part, the values of 5, and the modified geostrophlc wind 

components at z = /i (it^ and v/). To evaluate this dependence. 
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fields of both RAOB— and satellite-derived modified Ekman winds 
were computed using values of h ranging from 0 to 125 m AGL. 
Six values, at 25 m Intervals, were used In all. For /i = 0, a 
constant flux layer was not considered to exist, and Km. was set 
equal to Its classic value of 1.0 m^ s“^. Standard deviations 
of differences between "actual" winds and the RAOB and satel- 
lite Ekman fields were computed, along with corresponding 
linear correlation coefficients. Statistical results for wind 
speed are presented In Figs. 20 and 21. Recall that Gerrity's 
empirical model sets A equal to 50 m; this value Is supported 
by the RAOB data (Fig. 20) , yielding overall the smallest stan- 
dard deviations and highest correlations. Closer agreement 
between "actual" amd satellite-derived winds can be achieved by 
setting h equal to 75 n as indicated In the results displayed 
in Fig. 21. However, the 25 m increase in the depth of the 
surface layer produces "smoother" fields of satellite winds, 
which are more c<^patible with the observed data. This 
suggestion of a higher value for h is thus considered to stem 
frcan the same problem of differing degrees of detail outlined 
in the previous section. The value of 50 m, therefore, is 
deemed optimum for A. 

The second parameter under discussion is Km . Recall from 
Chapter 2 that is regarded as a function of height in the 
surface layer. Within the transition layer, however, the gvwan 
approach assumes the eddy viscosity to be equal to its value at 



Fig. 20 Vertical profiles of standard deviations (m s"*!) of the 
differences between "actual** wind speeds and RAOB-deri- 
ved modified Ekman wind speeds based on varying values 
of /i (left graph). Corresponding linear correlation 
coefficients are shown in the right graph. 
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z — h and Independent of altitude. In the current scheme, 

Is computed separately at each grldpoint In the plane z'’=h, 
Holton (1979) and Hess (1959) both note a mean value of '\-5.0 
s”^ forA",ni“ * stable boundary layer 1 km in depth. 

Holton's expression relating Km to the depth of the Ekman 
layer yields 13*0 m s for the current boundeu'y layer which 
extends to 1600 m AGL. Richardson numbers computed for the 
present case Indicate forced convective, though nearly stable, 
conditions exist throughout the surface layer; thus, a similar 
or slightly higher mean value for the turbulent exchange coef- 
f lolent should be anticipated ( see Appendix E for barocllnlc 
effects). Indeed, RAOB data yield a mean of 13.1 m^ s"^‘, while 
the average satellite-derived value of Km ia 13*6 m^ s“^. Hor- 
izontal depictions of the RAOB- and satellite-generated values 
of Km are shown in Figs. 22 and 23, respectively. Note that 
higher values (indicating greater instability) eu*e located in 
the western halves of both areas, while smallest values are in 
the vicinity of the eastern pressure ridge. 

Assigning an appropriate mean value to Km would simplify 
the computational procedure appreciably. A test of the slgnl- 
flcanoe of using individually computed turbulent exchange coef- 
ficients, rather than a mean value, was therefore indicated. 
Fields of modified Ekman winds were computed using a constant 
(mean) set, variously, at 1,0, 5.0, 10.0, and 15.0 m^ s~1 , 
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Fig. 23. Satellite-generated values of (m^ s“^). 
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Stand8u*d deviations of differences between these and "actual" 
fields were compared to previously derived statistics relating 
the "variable Ekman fields to the observed data. 

Corresponding linear correlation coefficients were also calcu- 
lated. Results of these comparisons, eigain represented by the 
statistics for wind speed, are shown In Figs. 24 and 25 for 
RAOB and satellite data, respectively. Comparisons Involving 
winds based on the "variable approach yield lower standard 
deviations and a higher degree of correlation. Of the constant 
values Investigated, 15.0 m s produced the closest agreement 
to fields derived fr«a a variable ATm* Although use of 15.0 
s~^ would provide comparable results at 600 m and 900 m AGL for 
both data sets, an even higher mean value would be needed to 
produce the variable /T^ pattern from 50-300 m AGL and above 900 
m AGL. Use of a (horizontally) variable Is thus indicated 
as the preferred technique. 

I 

Finally, the ability of the computational scheme to handle 
errors in satellite thermal data Is investigated. Recall from 
Chapter 4 (Section a) that the majority of past studies docu- 
menting these errors have presented comparison statistics which 
relate differences between corresponding RAOB and satellite 
soundings. A RMS error of 1.3®C Is representative of such stu- 
dies near the surface (Bruce it al. i 1977; Smith 1978; 

Schlatter, 1980; Scoggins gl.. , 1981). Also noted was a ran- 
dom error (or noise level) value of 0.5°C for satellite-derived 



Fig. Ik. Vertical profiles of standard deviations (m s~l) of the 
differences between "actual" wind speeds and RAOB-deri- 
ved modified Ekman wind speeds based on varying values 
of K„,(left graph). Corresponding linear correlation 
coefficients are shown In the right graph. 





Fig. 25. Same as Fig. 24, except using satellite-derived modified 
Ekman wind speeds. 
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temperatures alone, documented by Hillger and Vender Haar 
(1979) using structure function techniques, 

A test to determine the effects of these errors on the 
resulting wind fields was performed in the following manner. 
The original TIROS— N thermal data were perturbed using 
ccuaputer-generated random errors which were normally distri- 
buted about zero. F!ive sets of perturbed satellite tempera- 
tures were based on a standard deviation of 1.3®C and another 
five on 0,5°C, (Degree of cloud combination at a given sound- 
ing was not considered.) Each of the ten perturbed data sets 
exhibited a different combination of errors. As with the ori- 
ginal TIROS— N temperatures, these various data sets were then 
objectively analyzed using the Barnes (1973) scheme. Standard 
deviations of differences between the (objectively analyzed) 
original and perturbed thermal data are depicted in the upper 
half of Fig. 26, Note that in all cases, standard deviations 
of differences between the grldded data au*e less than 50 % of 
the deviations originally introduced at the sounding sites. 

Modified Ekman winds computed using the ten perturbed data 
sets then were compared to the unperturbed satellite-derived 
wind fields. Standard deviations of differences in wind speed 
are given in the lower half of Fig. 26. Winds below 600 m AGL 
were most affected by the deliberately perturbed temperature 
data, with greatest deviations near the surface. At 50 m AGL, 


fig- 26. Vertical profiles of standard deviations of the dif- 
ferences between original and perturbed TIROS-N ther- 
1^1 data ( C) (upper graphs) and their resulting mod- 
ified Ekman wind speeds (m s“l) (lower graphs). Stat- 
istics based on input random temperature errors of 1.3 
and 0.5®C are shown on the lefthand and righthand 

ffranhc: . T*0 GT>or* ^ -1 1 tr 
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a 1.3°C error in the Input temperature field yielded a differ- 
ence as great as 2.2 m s~^ from winds speeds based on unper- 
turbed satellite data. Assuming a 0.5°C error, on the other 
hand, resulted in maximum differences of less than 0.9 m s~^ at 
50 m AGL and less than 0.3 m above 600 m AGL. 

The sensitivity study indicates that the objective 
analysis scheme filters approximately half of the random error 
present in satellite temperature data. Of the two values 
tested, 0.5°C is considered to be more representative since it 
was deduced from satellite data alone. Thus, variations in 
wind speeds attributable to random temperature errors au?e 
expected to be less than 0.9 n s”^ near the surface and 0.3 m 
s"^ in the upper half of the PBL. 

e. Kinematic parameters 

Results of Sections b and c indicate that the "best" RAOB- 
and satellite-derived modified Ekman winds occur at 300 m AGL. 
Four kinematic parameters were computed from data at this 
level: horizontal moisture divergence, moisture advection, 
velocity divergence, and relative vorticlty. A fifth parame- 
ter, vertical motion, was determined at the top of the PBL 
using step-wise integration through the model levels below. 

Fields of horizontal moisture divergence were computed in 
three ways. First, a "control" data set was obtained using 
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RAOB-derlyed modified Ekman wlads and mixing ratios generated 
from RAOB dewpoint temperatures. A second set was computed 
using corresponding satellite-derived modified Ekman winds and 
mixing ratios. Finally, as a test of the TIROS-N moisture 
data, a third set was obtained by combining satellite Ekman 
winds with RAOB-derived mixing ratios. The first two pairs of 
fields are shown in Fig. 27. The combination of satellite 
winds with RAOB mixing ratios yielded patterns (not shown) 
nearly identical to those generated strictly from satellite 
data; differences were only ±0.1x10"® g kg"^ s"\ This indi- 
cates that, for this particular synoptic situation, horizontal 
gradients of satellite-derived mixing ratios are comparable to 
their RAOB counterparts in the lower portion of the PBL. 

The fields shown in Fig. 27 exhibit very similar features. 

Centers of moisture convergence (negative values) are located 

along the western and eastern borders of both areas, while 

moisture divergence is evident through the central section. 

Negri and Fonder Haar (1980) noted that moisture convergence 

>10"® g kg"^ s"^ was favorable for severe convective activity 

in their study utilizing mesoscale cloud- tracked winds. In the 

current case, however, maximum values of both RAOB- and 

satellite-derived moisture convergence are approximately 
—4 -1 • 1 

10 g kg s (western centers). Larger values were not 
expected as the data are synoptic-scale in nature and indicate, 
overall, weak dynamic forcing at this predawn period (Fig. 5). 



s *) generated from (a) RAOB- and (b) satellite-derived 
modified Ekman winds and moisture data at 300 m AGL. 
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Radar sunmarles Issued later in the day (Fig. 28) show 
that areas of moisture convergence agree closely with the 
intense convection which occurred in the region. Activity in 
northeastern Arkansas and the Missouri "bootheel" (Fig. 28, 
1135 to 1735 (31T) corresponds to centers of moisture conver- 
gence obtained from both sounding data sources (Fig. 27). 
Similarly, convective cells in southcentral Texas (already 
formed at 1135 GMT) occur in regions of convergence. Further, 
weather watch boxes Issued in the afternoon for the region 
extending from Nebraska to the Texas Panhandle correspond par- 
ticularly well to areas of both RAOB and satellite moisture 
convergence. It is noteworthy that the greatest moisture con- 
vergence based on satellite (as well as RAOB) data coincide 
with the strongest convective cells which developed in central 
Kansas by 2235 GMT. 

To further investigate the moisture divergence patterns, 
fields of its two components and ^z i^zw)) were com- 

puted. Values for Vz'i^iuj) (Fig. 29) were generally two orders 
of magnitude smaller than those for vu(% \^ 2 ) (not shown, but 
similar to Fig. 27) • Patterns of moisture advectlon at 300 m 
AGL show similar features for both satellite and RAOB data. 
Positive advectlon (recall sign convention of (5-8)) predom- 
inates in the western half of the area, with greatest values 
over Oklahoma ('''0.023x10"^ g kg“^ s“b . Negative moisture 
advectlon occurs eilong the extreme western and eastern edges of 



Fig. 28. Radar sununarles for 1135, 1435, 1735, 1935, 2035 and 
2235 GMT 19 April 1979. 






• 28 . Cont iriued . 




Fig. 28. Continued. 
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Fig. 29. Fields of moisture advection (10~^ g kg~^ s~^) at 300 
AGL derived from (a) RAOB and (b) satellite data. 
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the area, with an additional small center over the Kansas- 
Missouri border. 

Fields of horizontal velocity divergence were also com- 
puted (Fig. 30) . Since It contributes most of the variability 
in reflects the same general features exhibited 
by moisture divergence (Fig. 27)* Maximum convergence (nega- 
tive divergence) values of -3.3x10"^ s“^ and -2.9x10“^ s“^ are 
observed In western Kauisas for RAOB and satellite data, respec- 
tively. 

Patterns of relative vortlclty based on RAOB and satellite 
modified Ekman winds are given In Fig. 31 for 300 m AGL. 
Again, the two types of wind data produce very slmllcu* results. 
Both fields show positive vortlclty In the western portion of 
the grid 2 U'ea. Positive vortlclty is also indicated over 
southeastern Missouri and eastern Arkansas. These positive 
areas correspond to centers of moisture convergence detailed In 
Fig. 27. 

Finally, vertical motion was computed at the top of the 
PBL using the approach described In Chapter 5, Section g. 
Similar patterns result from the two data sets (Fig. 32). 
Specifically, maximum upward vertical motion (UVM) occurs In 
the western third of the grid area where the majority of the 
convective activity later developed (Fig. 28). Values derived 
from satellite data are slightly weaker than those based on 
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Fig. 30. Fields of horizontal velocity divergence (10 ^ s’**) at 
300 m AGL derived from (a) RAOB and (b) satellite data. 




Fig. 31. Fields of relative vorticity (10~5 s~l) at 300 m AGL 
derived from (a) RAOB and (b) satellite data. 






123 


Fig. 32. Fields of vertical motion (cm s“l) at the top of the PEL 
based on (a) RAOB and (b) satellite data. 
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RAOB data (e.g., 3»0 versus 3*6 cm s”\ respectively » in 
western Kansas). Weak UVM also is indicated by both types of 
winds in the southeastern portion of the network. Areas of UVM 
coincide closely with regions of moisture convergence (Fig. 
27) » velocity convergence (Fig. 30) and positive vorticity 
(Fig. 31) . This correspondence among parameters holds for both 
RAOB and satellite data. Finally, descending motion occurs in 
the central section where values are '''-1.0 cm s~^ in both 
oases. 

Again, it is prudent to recall that the satellite-derived 
vertical velocities emd other kinematic parameters au*e based on 
data observed approximately 2 h prior to their RAOB-derived 
counterparts. Given the already close correspondence between 
klnonatlc fields from RAOB and satellite sources*, it is con- 
ceivable that even better agreement would have been obtained 


from concurrent data. 


7. SUMMARY AND CONCLUSIONS 


In summary t study has examined the feasibility of 
using satellite^derived thermal data to obtain synoptic-scale 
winds in the planetary boundeu*y layer. A diagnostic computa- 
tional scheme based on the AFGWC-BLM modified Ekman equations 
was employed. Sixty-two predawn (0921 GMT 19 April 1979) 
TIROS-N soundings provided the data necessary to compute modi- 
fied Ekman winds at seven levels between the surface and 1600 m 
AGL. Similarly-derived wind fields based on 39 rawinsonde 
soundings (taken approximately 2 h later) were used as a basis 
for comparison. Comparisons were also made with "observed** 
winds generated from contact data for 26 of the RAOB soundings. 

Qualitative suid statistical comparisons revealed that 
satellite— derived modified Ekman winds show excellent aigreement 
with corresponding wind fields derived from RAOB thermal data. 
Satellite winds tended to be slightly weaker than their RAOB 
counterparts (due primarily to their slightly weaker thermal 
gradients), and exhibited a greater degree of cross- Isobar ic 
flow. Closest agreement between the two data sets was at 300 m 
AGL, the level at which nearly half of the RAOB soundings were 
affected by inversions. Although both sets of thermally-derived 
winds compared favorably with the "observed" patterns, the 
RAOB-derived values were slightly better in this reg 2 u?d. The 

Ekman winds were clearly superior to geostrophic winds in the 
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levels nearest the surface. 

A sensitivity analysis was performed to determine the 
effect on the computed winds of the depth of the surface layer 
(h) f horizontal variations of and the random error in the 
satellite-derived temperature data. Through comparisons with 
the "actual" wind fields^ h. s 50 m was found optimum for use in 
computing RAOB— and satellite-derived modified Ekman winds. 
Values of calculated separately at each gridpoint in the 
plane z - h f rather than the use of the mean value, also 
yielded more satisfactory results overall, particularly between 
the surface and 300 m AGL. And finally, variations in the wind 
speeds attributable to random temperature errors in satellite 
data were reduced (through objective analysis) to less than 0.9 
® near the surface and 0,3 m s“^ in the upper half of the 
PBL. 


Kinematic parameters computed using satellite thermal data 
were very comparable to their RAOB-derlved counterparts. Areas 
of horizontal moisture convergence, velocity convergence, auid 
positive vorticity were nearly coincident and aligned in 
regions which later developed significant convection. Moisture 
advectlon was shown by both satellite and RAOB data to be two 
orders of magnitude smaller than ved-ues of moisture 
divergence/convergence and, thus, relatively insignificant in 
this synoptic situation. Patterns of vertical motion at 1600 m 
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AGL that were deduced fr^ satellite winds in the PBL reflected 
values emd features similar to their RAOB-derlved counterparts. 
Regions of maximum upward motion (for both data sets) 
corresponded to the areas of greatest moisture convergence. 

The main conclusion to be drawn from this reseau'ch is that 
thermal gradients obtained from satellite data were sufficient 
to generate realistic synoptic-scale winds in the PBL. These 
wind patterns were comparable to those obtained from RAOB- 
derlved thermal gradients. The quality of the satellite- 
generated winds suggests that they could either augment, or be 
used in the lieu of (e.g. , over oceanic areas), RAOB winds in 
the lowest levels of the atmosphere. A further conclusion is 
that satellite-derived moisture fields, long considered to be 
poorly defined, yielded low-level patterns of horizontal mois- 
ture convergence which successfully indicated areas of later 

developing intense convection. 

Further research expanding upon the findings of this study 
are indicated in several areas. Two limitations have been the 
smallness of the size of the sample case and the fact that only 
one map time was available. Generation of satellite-derived 
modified Ekman winds (and associated kinematic parameters) at 
several time periods would provide a temporal continuity not 
available during the current study. VAS satellite soundings, 
taken at 3 h intervals for example, would be particularly 
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suited to such research. Additional case studies, Including 
the selection of a more unstable period, would further test the 
usefulness of the computational procedure. 

The aim of the present study was not to explore the intri- 
cacies of the AFGWC-BLM (e.g., systematic model errors) but 
rather to use this operational model as a tool to medce prelim- 
inary calculations of satellite-derived PBL winds. It is felt, 
however, that the model could yield improved results by allow- 
ing the eddy viscosity to vary with altitude, e.g., linearly 
decreasing /T,, with height such that it vanishes at the top of 
the Ekman layer (as suggested by Estoque, 1973). Application of 
satellite data to other (improved) boundary layer models which 
deal with mean motion is also worth future examination. 


APPENDIX A 


DERIVATION OP MODEL EQUATIONS 
DERIVATIONS OF SURFACE LATER EQUATIONS 

Recall equations (3-3) f (3-4), and (2-11): 








11 + ^ 


-2 



kz 

1-fi 

J. 

ds 





if 

dz Cp 

3z 





Urn 


Taking the square root of (A-7): 


K^kz 


1-/S* 


ar 


dz 


ds 

dz 


-2 


Urn 


= kzu* - kziu 4S- 

11 + ^' 

ds 

d 

dz Cp 

dz 


-2 


Rearranging: 




ds . ds 
~ = itzit* — 

dz dz 


■0 


dT 

dz 


Hultl plying again by Ki 

18 


1 8 f H 

~ dz 


fczti* - fczii» 4?-Aj4^ + 


dz 


Substituting (A-1) and (A-2) Into (A-10); 

Ids 


ru* = XL? 


dz 


kzu» — kzu»-^M-ut-d* 
1? 


Dividing by uh 


u? = 


ds 

dz 


\kzu» — fcziJ. 


U 

■d 


Rearranging: 


ds 


i:zu* 


= tx.* + kz-d,^ 
if 


Dividing by kzxu: 


ds _ J?2. 

dz ~ fcz u. >iJ 
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(A-7) 


(A-8) 

(A-9) 

(A-10) 

(A-11) 

(A-12) 

(A-13) 

(A-14) 


Integrating with respect to z from z, to ht 


(A-15) 


^ y 1 V Ba 

= — /* (-)dz + — — /I 
k'^*» z It* 


s(z-h) - s(z-z^) = ln(-^) + ^ -^(/i - 2o) 

K Zq U* ^ 

The wind epeed Is assumed to vanish at 2 = 2 ,, thus: 
S - .(»=*) = ln(A) + £1 m.(h - z.) 


HylYattgB 2 .: Forced 
Substituting (A-1) into (A-2): 




ds 

-1 

11 + ^ 

dz 


dz Cp 


Rearranging: 



ds 

dz c. 

3z 


-1 




— 

u» 


Substituting (A-1 8) into (A-3); 

1-1 


«.=! 


ds 


dz 


— 

Urn 


Rearranging (A-1): 




ds 


dz 


-1 


Substituting (A-20) into (A-19): 

^ 5 u/ u. 

Inverting (A-20) and substituting into (A-4) 
JC^[kz(l-fiR,)f^ 
Multiplying (A-22) by 

XL? 


(A-16) 


(A-17) 


(A-18) 


(A-19) 

(A-20) 


(A-21) 


(A-22) 
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^ s [kz{l-pRi)f (A-23) 

tx* 

Taking the square root of (A>23): 

^ = (A-24) 

Substituting (A-24) Into (A-21): 

Ri=^^ kz(l-pRi) (A-25) 

V XL* 

Let us now ecamine the term (l-fiRi): 

(l-fiK,) . 1 - (A-28) 

Substituting (A-25) for the last Ri on the RHS of (A-26): 


« 1 - (l-p/?i) ^^kz (A-27) 

tJ XL* 

Rearranging . ( A-27 ) : 

(,1-pRt) + ^^kz = 1 (A-28) 

(1-^) [1 + ^ ^ fcz ] = 1 (A-29) 

(1-^) = [1 + ^ fcz ]-^ (A-30) 

u** 

Substituting (A-30) into (A-24): 

if = kTXL*[\ + 42- fcz]“' (A-31) 

1? u» 


Evaluating K zt z = h, the exchange coefficient for momentum 
for forced convection conditions becomes: 

= fc/Mx.[l + 42- ^ jkA]-‘ (A-32) 

1 ? u/ 

[As an Intermediate result of the above derivation, 
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(A-30) may be rearranged to yield the alternate expression 
for Ri given In equation (3-19).] 


£Cfi& Convection Regime 


Recall the "free" expression for the eddy exchange coeffi- 
cient, equation (3-5): 

^ I ^ ^ (A-33) 

where X Is a dimensionless constant equal to 1,2 (Gerrlty, 
1967). In terms of momentum, however, the "free" eddy viscos- 
ity Is only 70% of (A-33) (Priestly, 1959), l.e., 

(A-34) 


0.70 

! 


91 + JL 

1/^ 


1 

1? 

dz Cp 



Derivation 2: Free ^ 


Rearranging (A-2), take the absolute value (recalling that A 
Is always positive): 


+ JL 
dz c, 


K 


Substituting (A-35) Into (A-34): 


JC - 0.70 


|u.-d. 1 

1/2 



K 



Squaring (A-36): 


X* - 0.49 s* i 

•tf AT 


(A-35) 


(A-36) 


(A-37) 


Multiplying by Kt 


K* = Q.A9 z* 


(A-38) 
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Rearranging (A-1) and cubing; 


Substituting (A-39) into (A-38): 


■u»® = 0.49 X® z * ^ I u.iJ. 


Rearranging: 


*2.04 


Tadcing the cube root: 


ti»* -a 


X* z* g |u»‘d»| 


0* r x*z* y 


* (2.04)*'^ 


= 1.27 


« r-it/3 




r—ii/s 


Gerrity notes that between z, and (zo + lm.) a logarithmic 
wind profile 1s assumed to exist. By definition, the loga- 
rithmic wind profile for turbulent flow is given by: 


u, z+z„ 

s = In 

k z. 


(A-43) 


where z and z, are comparable; for our purposes, z = l m. 


Thus, 


/■•*' is-d. - T = in ^ 
•'*. a* fc 2 - 


. . , t «o + 1 

S(zsz«-K) - s(zsz«) = -7T 1" — 


(A-44) 


Now integrating (A-42) with respect to z from Zo + lm. to hi 


135 


A '« = 

S 1-27 


u* 




\u,i}.\^'^ 

1/18 


1/3 




dz 

/.^l 


S 1.27 


tx/ 


X«^ Iwtf, 1.*'^ 

s(zs:X) - s(z=Zo + l) s -3.8 


y 


1/3 


ti,» 


[-3z-‘^]L., 


— ]l/» 

T> 

X 

91 


X«^ |u,T>,|*^ 

[A->/3-(2,+1^)-./3J 

Substituting (A»44) into (A-45) and rearranging: 


(A-45) 


S..(*=A) = :^l„^]-3.8 


tu^ 




1/3 


(A-46) 


•Derivation 4: Free 


An expression for "free" K„ is obtained by following the 
above derivation of .? to equation (A-38). Evaluating (A-38) 
at z = A and taking the cube root: 


0.49 x*a“ 

i 

1 ix«t9. 1 


1> 



The constant is conmon to both free and forced convec- 
tion equations. Deriving an expression for using the "free" 
equations yields a highly constrained equation involving com- 
plex variables. However, a more suitable expression may be 
obtained using the "forced" equations. 

Derivation «:♦ 


Hearrsngiog ( A— 1 ) 2 


136 


ds 

dz 


lU 

~K 


Rearranglog (A-2): 

Substituting (A-49) Into (A-48): 


dz Cp 


-i 


ds 

u#* 


+ -2- 

dz 

” tu'a. 

dz 

*r 

J 


Recall (A-1A). Substituting (A-50) into (A-14): 

M*[ar 


+ -2- 
iU[dz Cp 


-}^ + 2Lia. 

kz u* -a 


(A-48) 


(A-49) 


(A-50) 


(A-51) 


Multiplying by — : 

XU 


1. 
u 

Qz Cp 


-±L + i!.i£ 

~kz " “ 


u/ 1? 


(A-52) 

Integrating with respect to z from the height of the instrument 
shelter, to A: 


+ jz_ 

dz - L 

iL^^Lig. 

dz ^ 


kz u? 1 ? 


dz 


r -^<£* + r -^dar = f -r-dz + f — g- ^dz 

dz •'*» Cp •'■i fcz •'•i It. a 


(Th ~ Ti) + ~ ~fc~ ~ ^ (*~2i) 




u*^ a 


(A-53) 

(A-54) 

(A-55) 


Dividing by (A-z<) and rearranging; 


«tt.” 


to(^) 

*i 






(A-56) 


A Zf >,p 

Equation (A-56) can be solved for i>. through use of the qua- 
dratic formula. Let: 


a^JS- 

•&U* 


b = 


fc(A-*i) • 


c = — 




h—Zi 
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Thus, 

+ 6i?« + c = 0 

Using the quadratic foniiula: 

- _ —6 ± v6®-4oc 


There 2 u*e two possible roots to (A-58): 


—6 + Vb®— 4oc 


and r. = r& ~ 

* 2o 


i9« must vanish for neutrea stratification. Under neutral con- 
ditions the lapse rate becomes — SL ; looking at the constant 




dz Cp 


+ s_ il + JZ. = 2_+ ^ = 0 


Thus, Ti Is the appropriate solution, giving: 

1 I r .. tr r 1 

fi_ + + _2_ 

m 


2flg k{h-Zi) I*(h-Z 4 )j ^2 K-Zi Cp 


2fg *(*-.,) 


,/g ■°<i> 

2fig k(h-Zi) 


k(h-Zi) 

ln(^) 

lo(^) 

*i 

*(*“*<) 

1.0+iSi 







(A-59) 


TRANSITION LAYER EQUATIONS 


PerlYatlOB 1: Modified B»nnan Equations 


Recall froB Chapter 3 (Sections a and d), the horizontal 
equations of motion, (3-1) and (3-2), and the simplifying 
assumptions (3) and (4): 




12. +oa— 

+ _a_ 

K — 


p 

dx dx 

dz 

^"*azl 


0 * /v - fVg + 




and 


dv ^ l|dp . __9^ 
-TT- = -/t* + pq-r— 


dt 


P |9v 


pg 




dy 

dz' 


d [y av ] 
dz dz\ 


(A-60a) 

(A-60b) 

(A-61a) 

(A-61b) 


He now Introduce the complex variable i , where i = and 
=: -1. Multiplying (A-6lb) by i and combining with (A-60b): 


0 = -i/u + t/Ug + Kn (w ) (A-62) 

02 

0 = fv - fVg + M 

■ p 

0 = -ifu + fv + ifUg - fVg + Kn — Y (u + iv ) (A-63) 

dz 

Rearranging: 

a2 

0 a i/(ii, + Wg) -i/(u + w) + + w) (A-64) 


Now we define the complex variables: 

V a u + tv and Vg = Ug + ivg 
Substituting these Into (A>64): 

a^v 
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W ss o + 6 (A-70) 
We now have two boundary conditions to consider: the upper and 
lower eztreses of the transition layer. Classically, these 
limits are set equal to Infinity («») and zero, respectively. 
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First : When then V-»V^ and W-*0, Setting ff 

equal to <» and plugging these values Into (A»70): 

0=0 et***^" + b e ■!*•*'*)*"“ 

Since e“ = " , this .Implies that a = 0 and we are left 
with: 

W = 6 (A-71) 

Second B.C. : When z-»A, then W(=V-V,) -♦ where V Is 

the wind computed from applicable surface layer equations 
(9 s [/ + iy) ^ and Is the geostrophio wind at z = h 
(V^ = + ivg). To set the height coordinate equal to zero 

at Af let us redefine z as (z-h). Plugging these values 
Into (A-71): 

^-V^= 6 g-«+l)mO 

Since e° = 1, 6 = V-V^ and the particular solution 

becomes: 

W = (V-Y^) e (A-72) 

Expanding the various terms, (A-72) becomes: 

V-Va = (V-V^) 

(u*iv) - (ii,+ivg) = [(£7+ iV) -(i4 + iv^)] (A-73) 

Applying Euler's rule (e^^ = cos d + isin i?) to (A-73) : 

(u+w) “• (ltf+ iVg) 

= i(U + iV) - {tig + tWj)] e"”“ [oos(mz) - i sin(mz )] 

= (£7+ iV) [cos(t7iz) — isin(mz)] 

- (vg + Wg) e~^ [cos(mz ) - isin(mz )] 
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= Um~^ [cos(mz) - isin(mz)]’ 

+ iVe~^ [cos(mz) — isin(mz)] 

— tig [cos(ma: ) — isin(mz )] 

— ivg e""**[cos(mz) - isin(mz)] 

= Ue~^coa(mz) - iUe~^^n(mz) 

+ iVe~*^ oos(7m) + P'e~~’sin(mz) 

- «g e"*** cos(t7iz ) + ii£^ e~"** sin(7nz ) 

- «"~coa(m*) - e"*~sin(mz ) 

(A-74) 

Separating the real and imaginary . parts of (A-74), we obtain 
the following equations: 

u-ti, a (7«"~coa(mz) + Ve~^sm(mz) 

- ti* coa(mz ) - e“"“sin(mz ) 

- • "** eoaimz) + (V-Vg) sin(mz)] (A-75) 

v—Vg ~ —</€""•* ain(mz) + J'c~"“cos(mz) 

+ Ug«“*“sin(mz) - Vg e””*cos(mz) 
a e"^ [(V-t»^) coa(mz} - (U-u^) sin(mz)] (A-76) 

Finally, recall that z was redefined as (z-h). Making the 
appropriate substitutions and rearranging: 

u a tf, + cos[m(z-h)] + (7-v^) sin[Tn(z-/i)1J 

/ . (A-77) 

V a V|g + •”"* * t/g)co3[m(z— h)] — sin[m(z— h)]j 

(A-78) 

OarlTattOB X: Modified Geoatroohle Eouatlona 

Recall assumption (5) from Chapter 3, Section d; as noted 
previously, if is linear with height, then the change of 
geostrophlc shear with height must be zero: 
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The general solution to (A-79) Is of the form: 

= o + 6z (A-80) 

Before determining a and b, let us redefine z as (^-z) where 
h'^ z ^ H. Equation (A-80) then becomes: 

V, ==a + b{H-z) (A-81) 

We now have two boundary conditions to consider: 

First B.C. : When z-*H^ then Plugging these into 

(A-81): 

+ b(H-H) 

Thus, a ^ and (A-81) becomes: 

^g ^P“*b{H-z) (A-82) 

Second B.C^ : When z-»A, then Plugging these into 

(A-82): 


6(/f-A) 

Thus, 6 = ~ and the 


solution 


becomes: 




This equation is not in a form that suits our needs. We want 
the component form of (A-83) expressed in terms of temperature. 
To obtain a more usefiU. form of (A-83), 1st us consider first 
the Hypsometric Equation: 


*1 = In 

9 

RT^ , 

= In 


P2 


Pi 


Pi 


^2 



where pz< P] 
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Letting 1 = /i and 2 = H, and rearranging; 

-=1 g = ^ 

r* j {ph j 

Differentiating (A>81>) with respect to x; 

te [I-fI * ''"'H " 1 


In 

Pa 


^ , Pa 
+ In 


Pi^ 


PH 


_ 

I Pa [Ph 



1 

9pA 

Ptf 

dx 


1 ^Ph 9pff 

Pa 9* Ph 9a: 


From the Equation of State we can write; = p^RT^ and 
Ph=Ph^Tu. Substituting these into (A-85); 


d 

1 

17 (*ff-*A) 


1 

9PA Ptf 

1 

^Ph 

dx 


“ PaT’a 


dx pu Tu 

Ph 

dx 



1 

dpH 

PhTh 

dx 


Multiplying (A-86) by 7*^ 

I U-1 »<*--. 




1 dPH 

Ph Tfr 9x 


Recall that: Vg = ^ ^ Assuming constant density: 


-= ir (*it-*a) = [/«* - 


dE 

Th 

dx 

Th 


Dividing by / and rearranging: 
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’■* / to [ * / to h^l * 

Subtracting tvcm both sides and rearranging terms (b), (c) 


and (d): 


(c) (d) 


te * / to < 

... I g ag[rj] g aglri.1 

f dx Tu f dx Ta " T;/ 


- V? = - V? 


.H 


g dE\Th Tg a Zh_ - IlL 

/ axfTV " 7WJ ^“[Tg Tg\- 

... .JL I + = 


/ ar Tg 


Rearranging: 


r^iL -g,[ [-^l/g )[ - 1. I .yullltZllL + 

/ axllr*)''^ *'] / ax[ Tg ) ) 


Differentiating (A-84) with respect to y: 


,Jf -..h 


(A-91) 




(as in (A-85)) 

I fl(j. 1 (A-9; 

Ipj* ay ph ay J 

Following the same steps as in (A-86) through (A-91), and 
recalling that tig = — ^ ^ ^ ^ vre obtain: 

p/ ay / ay 


-n-2- J- 4= (Zg-Z„) + ^ 


dE T/^-Tg 

ay T’y 


Th 


+ u^ = 


(A-93) 

Returning now to our particular solution, (A-83) can be broken 
down into component form: 


Solution £sz u^: 
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Substituting (A-93) ^or and recalling that Zh - 

Zif = Hi 




Th-Th 


Th 


a 

1 1 

ay 

T*l 


• ~ (7 dE 

where: u = — — 

/ dy 


SsliiiiSR £sl 


Vm — V„ + — — — (Vn ~ Va) 

Substituting (A-91) for v^\ 

\Tk-Th 


V, = wf + 


!<vf - v) 






r* 


and 

(A-94) 

(A-95) 

(A-96) 


where 


- o dE 

V = -“ 

/ dx 


(A-97) 


APPENDIX B 


RESPONSE DETERMINATION 

The Barnes (1973) objective analysis scheme requires 
selection of two weight function constants (g and 4c) to 
achieve the desired eunplltude response at specific wavelengths. 
Input temperature, dewpoint, and pressure data were analyzed 
using the ten c<^binatlons of g and 4c given in Table 10. The 
respective response curves have been shown in Fig. 8. The 
resulting RAOB- and satellite-derived modified Ekman winds were 
then compared statistically to ’’actual" wind fields (objec- 
tively analyzed at corresponding responses) to ascertain the 
most appropriate degree of resolution. The objective analysis 
that yielded the lowest standard deviations of differences 
between "actual” and modified Ekman winds, as well as the 
highest linear correlation coefficients, was Judged most 
appropriate for the investigation. Results from three analyses 
('^50} response at wavelengths of 900, 1000, and 1100 km) which 
most closely matched these criteria are shown in Figs. 33 and 
34. Of the three, 900 km exhibited the "best" overall results 
using RAOB data, but the "worst" results when satellite winds 
were employed. Conversely, 1100 km had the opposite effect. 
Thus, a response yielding 50{ resolution of amplitude at 1000 
km was deemed best for the three types of data ("actual", RAOB, 
satellite) . 
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Table 10. Weight function constants selected to achieve approx- 
imately 50% response at the indicated wavelengths 
using the Barnes (1973) objective analysis scheme. 



4c 

Wavelength (km) 

0.2 

30,000 

300 

0.4 

50,000 

500 

0.3 

120,000 

700 

0.4 

160,000 

900 

0.4 

200,000 

1000 

0.4 

240,000 

1100 

0.4 

400,000 

1400 

0.4 

450,000 

1500 

0.4 

575,000 

1700 

0.4 

800,000 

2000 
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differences between “actual" wind speeds and RAOB-de- 
rived modified Ekman wind speeds objectively analyzed 
for a 50% response at wavelengths of 900, 1000 and 1100 
km (left graph) • Corresponding linear correlation co- 
efficients are shown in the right graph. 
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APPENDIX C 


CORRELATION COEFFICIENTS 
LINEAR CORRELATION COEFFICIENT 


The product-moment formula for the linear correlation 
coefficient, r, is (see, for example, Panofsky and Brier, 


1968 ): 

2 _ 

(Ex2)(Ey2) 


(A-98) 


where x = Jf-P’andy = Y-Y. AT and K represent the i-th value of 
two scalar quantities (temperature, wind speed, etc.) being 
compared. 


VECTOR CORRELATION COEFFICIENT 

Court’s total vector correlation coefficient, /?, relates 
the wind at one time or space to wind at another time or space 
(Lambeth, 1966). Given two samples of wind: 

^ = Ux + Vj and 2 = Xx + Yj 

the vector correlation coefficient is obtained through the 
expression 

(s2+Sv")(s.V-4) 

Here u, v, x and y are deviations from the means, i.e., 
u = U - U, etc. Tlie terms s*, and denote sample 

variances of u, v, x and. y,^ respectively; s„*, s„», Md 

Sxy eu*e covcU?iances of the indicated deviations. 
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APPENDIX D 


5CALK AMALYSIS DEL HORIZONTAL MOMENTUM EQUATIONS 

The appropriateness of neglecting individual derivatives 

(and, hence, the advectlve terms) in the horizontal manentum 

equations (Chapter 3» Section d) can be seen in the scale 

analysis given below. Since only one map time was available 

for the current study, and were assumed to equal zero. 

ot at 

Values of the various terms were computed using "observed" data 
at 50, 300 and 1200 m AGL for six points representative of the 
research area (Fig. 35). Averages of values at the six points 
produced the results shown below. Horizontal and vertical 
advectlon terms were one to two orders of magnitude smaller 
than the Coriolis, pressure gradient and frictional terms at 
all three levels. 



fv, — /u: ~10“^ m s"* 


-ilE. -10-3 „ 

p ax p ay 


, „ d^v 

"‘as*- 


m s"* 
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Fig. 35. Data points used in scale analysis. 


152 



APPENDIX E 


BAROCLINICITY IN THE PEL 

The PBL under study is noted In Chapter 6, Section d, as 
being forced convective in nature, though nearly stable. The 
weak barocllnlclty that is implied can be described by the 
amount of vertical geostrophlc shear and, near the surface, the 
degree of stress and cross-isobarlc flow. 

The terms geostrophlc shear and barocllnlclty may be con- 
sidered synonymous. Inserting RAOB-derived thermal data into 
(3-3^) and (3-35) (forms of the thermal wind equations), a mean 
geostrophlc wind shear of 0.8 m s"^ km"^ is obtained over the 
reseeurch area. This is much smaller than is the ''-2.0 m 
km"^ shear given by Krishna (1981) and is less than one-fourth 
the 3.5 m s"^ km~^ noted by Arya and Wyngaard (1975) for a typ- 
ical midlatitude barocllnic PBL. 

The amount of stress ( r) in the surface layer may also 
indicate the degree of bau*oclinicity. Values of surface stress 
tabulated by Hoxit (197*1) for an average baroclinic PBL range 
froD 0.56 to 1.05 N m“^. In the current study, however, values 
at the top of the surface layer are much smaller, with a mean 
value of '^0.13 M m”^ for each of the RAOB- and satellite- 
derived modified Ekman and "actual" winds. Since stress in the 
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surface layer is considered to be constant with height (Hess, 
1959; Holton, 1979) » the study values are assumed to be 
representative of stress at the surface. Stress profiles above 
the surface layer are shown in Fig. 36 for the three wind data 
sets. (t is obtained through (2-1).) "Actual" stress values 
indicate that the region of maximum barocllnicity occurs at 300 
m AGL, the average height (top) of the low level inversion. 
This maximum (0.157 N m"^), however, is also significantly less 
than Hoxlt's values. The RAOB- and satellite-derived profiles 
show a constant decrease with height, as is "ideally" expected 
to occur (Hess, 1959). 

The degree of barocllnicity can also be inferred from 
cross-isobaric flow. Ekman winds (as discussed in Chapter 2, 
Section c) flow toward low pressure at or near the surface, 
making an angle of 45° with the isobars. The modified Ekman 
equations, however, partially allow for the effects of baro- 
clinicity (via the temperature gradient terms in (3-34) and 
(3“35)). Thus, a decrease in the cross-isobaric angle at or 
near the surface is anticipated. In the current study, mean 
cross-isobaric angles at the top of the surface layer (50 m 
AGL) 8u?e 26.9° (with air flowing toward lower pressure) for 
RAOB-derived modified Ekman winds, and 28.6° for their 
satellite-derived counterparts. "Actual" winds near the sur- 
face yield a mean cross-isobaric flow angle of 33.9°. Strong 
barocllnicity should yield much smaller crossing angles as the 



Fig. 36. 


Stress profiles based on RAOB- (solid) and satellite- 
derived (dashed) modified Ekman data, as well as 
actual" (dotted) winds. 
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slurface winds would be forced (through instability) to flow 
more parallel to the pressure contours. 
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